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Given a class F{6) of differential equations with arbitrary element 8, the problems of sym- 
metry group, nonclassical symmetry and conservation law classifications are to determine for 
each member / G J~(8) the structure of its Lie symmetry group Gf, conditional symmetry 
Qf and conservation law CLf under some proper equivalence transformations groups. 

In this paper, an extensive investigation of these three aspects is carried out for the 
class of variable coefficient (l+l)-dimensional nonlinear telegraph equations with coefficients 
depending on the space variable f[x)uu — (g{x)H(u)u x ) x + h(x)K(u)u x . The usual equiva- 
lence group and the extended one including transformations which are nonlocal with respect 
to arbitrary elements are first constructed. Then using the technique of variable gauges of 
arbitrary elements under equivalence transformations, we restrict ourselves to the symmetry 
group classifications for the equations with two different gauges g = 1 and g = h. In order 
to get the ultimate classification, the method of furcate split is also used and consequently a 
number of new interesting nonlinear invariant models which have non-trivial invariance al- 
gebra are obtained. As an application, exact solutions for some equations which are singled 
out from the classification results are constructed by the classical method of Lie reduction. 

The classification of nonclassical symmetries for the classes of differential equations with 
gauge g — 1 is discussed within the framework of singular reduction operator. This enabled 
to obtain some exact solutions of the nonlinear telegraph equation which are invariant under 
certain conditional symmetries. 

Using the direct method, we also carry out two classifications of local conservation 
laws up to equivalence relations generated by both usual and extended equivalence groups. 
Equivalence with respect to these groups and correct choice of gauge coefficients of 
equations play the major role for simple and clear formulation of the final results. 

Mathematics Subject Classifications (2000): 35LI0, 35A22, 35A30 

Keywords: symmetry classification, nonclassical symmetry, conservation law, equiv- 
alence group, nonlinear telegraph equation, exact solutions, symmetry analysis, Lie 
algebras 



1 Introduction 

Since the notation of continuous group was introduced by Lie at the end of 19th century, signif- 
icant progress in application of symmetries to analysis of concrete nonlinear differential equa- 
tions has been achieved. The classical Lie symmetries of nonlinear differential equations allow 
us to find explicit solutions, conservation laws, linearizing substitutions of the Hopf-Cole type, 
etc [I3 l l43 p44p70p73| . For most of the application the search of explicit symmetry structure of the 
corresponding differential equations is always a crucial step, which consist of the cornerstone 
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of group analysis of differential equation. Consequently, there are two interrelated problems 
which still remain to be solved in the traditionally group analysis of differential equations. The 
first problem consists of finding the maximal Lie symmetry group admitted by a given equa- 
tion. The second problem is one of classifying differential equations that admit a prescribed 
symmetry group. The principal tool for handling both problems is the classical infinitesimal Lie 
method [13, 43, 70, 73j. It reduces the problem to finding the corresponding Lie symmetry alge- 
bra of infinitesimal operators whose coefficients are found as solutions of some over-determined 
system of linear partial differential equations (PDEs). However, if the equations under study 
contains with arbitrary element (functions or parameters) , then one has to solve an intermediate 
classification problem. Namely, it is necessary to describe all the possible forms of the functions 
involved such that this equation admits a nontrivial invariance group. Generally, the problem 
can be described as follows: Given a class J-{9) of differential equations with arbitrary element 
9, the problem of group classification is to determine for each member / G F(9) the structure of 
its Lie symmetry group Gf, or equivalently of its Lie symmetry algebra Af under some proper 
equivalence transformation groups. This description is also fit for the problems of nonclassical 
symmetry and conservation law classification by replacing Lie symmetry with these two different 
notations. 

Historically, the first classification of Lie symmetries was derived by Lie, he proved that a 
linear two-dimensional second-order PDE may admit at most a six-parameter invariance group 
(apart from the trivial infinite parameter symmetry group, which is due to linearity) [62J. The 
modern formulation of the problem of group classification of PDEs was suggested by Ovsiannikov 
in 1959 |72j . in which he present complete group classification of a class of nonlinear heat 
conductivity equations by using the technique of equivalence group and direct integration. After 
that, the group classification of nonlinear PDEs became the subject of intensive research. A 
detailed survey of the work done in this area up to the beginning of the 1990's is given in |44j . 

In the present paper we investigate a class of hyperbolic type variable coefficient (1+1)- 
dimensional nonlinear telegraph equations of the form 

f{x)u tt = {g{x)H(u)u x ) x + h(x)K(u)u x (1) 

where / = f(x), g = g(x), h = h(x), H = H(u) and K = K(u) are arbitrary and sufficient 
smooth real-valued functions of their corresponding variables, f(x)g(x)H(u) ^ 0. In what 
follows, we assume that (H U ,K U ) ^ (0,0), i.e., ((TJ is a nonlinear equation. This is because 
the linear case of (JI]) (H, K = const) was studied by Lie |62| in his classification of linear 
second-order PDEs with two variables. (See also a modern treatment of this subject in |73j). 

The study of equation (pQ) is strongly stimulated not only by its intrinsic theoretical interest 
but also by its significant applications in Mathematics and Engineering. In fact, hyperbolic type 
second-order nonlinear PDEs in two independent variables are usually used to describe different 
types of wave propagation. They are also used in differential geometry, in various fields of hydro- 
and gas dynamics, chemical technology, super conductivity, crystal dislocation to mention only 
a few applications areas. The corresponding models are comprised by the Liouville, sine/sinh- 
Gordon, Goursat, d'Alembert, Tzitzeica and nonlinear telegraph equations and a couple of 
others. From the group-theoretical viewpoint the popularity of these models is due to the fact 
that they have non-trivial Lie or Lie-Backlund symmetry [IOllTgll2lfl l£fl ^l7Tni751IH5] . By this 
very reason they are either integrable by the inverse problem methods or are linearizable and 
completely integrable [T|l2"6l 66J. 

The investigation of Lie symmetry classification of the (1 + l)-dimensional hyperbolic type 
second-order nonlinear PDEs has a long history. Probably, Barone et al [8] was the first study 
of the following nonlinear wave equation u%t = u xx + F(u), by means of symmetry method, 
this equation was also studied by Kumei [56] and Pucci et al [83] subsequently. Motivated by 
a number of physical problems, Ames et al [HE] investigated group properties of quasi-linear 
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hyperbolic equations of the form 



u tt = [f(u)u x ] x . 



(2) 



Later, their investigation was generalized in |23 [ l46 t l57] to equations of the following forms 
respectively 

u tt = [f(x, u)u x ] x , u tt = [f(u)u x + g(x, u)] x , and u tt = f(x, u x )u xx + g(x, u x ). 

The alternative form of equation ([2]) was also investigated by Oron and Rosenau |71j and Suhubi 
and Bakkaloglu [86]. Arrigo [7j classified the equations uu = u™u xx + F(u). Furthermore, 
classification results for the equation u u + K{u)ut = [F(u)u x ] x can be found in |l4"tl71|. An 
expand form of the latter equation uu + K(u)ut = [F(u)u x ] x + H(u)u x was studied by Kingston 
and Sophocleous [55] , Recently, Zhdanov and Lahno [59] presented the most extensive list of 
symmetries of the equations 

uu = u xx + F(t,x,u,u x ) 

by using the infinitesimal Lie method, the technique of equivalence transformations and the 
theory of classification of abstract low-dimensional Lie algebras. There are also some papers [W\ 
[32j[7H[82] devoted to the group classification of the equation of the following form 

u tt = F(u xx ), u tt = F(u x )u xx + H(u x ), and u tt + Xu xx = g(u, u x ). 

It is worthwhile mentioned that the constant-coefficient nonlinear telegraph equations 

u tt = {F{u)u x ) x + H{u)u x 

together with its equivalent potential systems have also been studied by Bluman et al [Iipi4ffl6] . 
In their a series of papers, many interesting results (especially for case of power nonlinearities) 
including Lie point and nonlocal symmetries classification and conservation law of the four 
equivalent systems were systematically investigated. Recently, Huang and Ivanova present a 
strong complete group classification for a class of variable coefficient (l+l)-dimensional nonlinear 
telegraph equations of the form [37] 



Exact solutions and classifications of conservation law with characteristics of order were also 
investigated [37] . 

From the above introduction, we can see that equation (pQ) is different from any aforemen- 
tioned ones and is a generalization of many well studied equations. What's more, equations ([I]) 
can be used to model a wide variety of phenomena in physics, chemistry, mathematical biology 
etc (see Section 2 for detail) . Thus there is essential interest in investigating them from a unified 
and group theoretical viewpoint. 

In this paper, extended group analysis of class (1) is first carried out. The usual equivalence 
group and the extended one including transformations which are nonlocal with respect to arbi- 
trary elements are constructed for class (1) and its subclasses. The structure of the extended 
equivalence group and non-trivial subgroup of (nonlocal) gauge equivalence transformations are 
investigated. As a result, group classification problems related to two different gauges g = 1 
and g = h are really solved for each form with respect to the corresponding usual and extended 
equivalence group. Classical Lie reduction of some classification models are described and exact 
solutions are obtained by using the reduction. Nonclassical symmetries classification of class (1) 
with the gauges g = 1 is discussed within the framework of singular reduction operator. Several 
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nonclassical symmetries for equation form class (1) are constructed. This enabled to obtain some 
exact solutions of the nonlinear telegraph equation which are invariant under certain conditional 
symmetries. Using the most direct method, two classifications of local conservation laws up to 
equivalence relations are generated by both usual and extended equivalence groups. Equivalence 
with respect to these groups and correct choice of gauge coefficients of equations play the major 
role for simple and clear formulation of the final results. 

Problems of group classification, except for really trivial cases, are very difficult. Generally, 
there are two main approaches in studying group classification problems in the literature. The 
first one is the algebraic methods and is based on subgroup analysis of the equivalence group 
associated with a class of differential equations under consideration. Its main idea rely on the 
description of inequivalent realizations of Lie algebras in certain set of vector fields of the equation 
under consideration [9,93], which was original from S. Lie [44,62] and recently rediscovered by 
Winternitz and Zhdanov et al [34^93], The method has been applied to classifying a number of 
nonlinear differential equations [2 ,'9, 33, 34, 40-42 ,59 ,60, 93-95], including the class is normalized 
(see [81] for rigorous definitions of normalized classes and related notions). The second approach 
is based on the investigation of compatibility and the direct integration, up to the equivalence 
relation generated by the corresponding equivalence group, of determining equations implied by 
the infinitesimal invariance criterion [73]. This method was suggested by L.V Ovsiannikov and 
referred as the Lie-Ovsiannikov method. This is the most applicable approach but it is efficient 
only for classes of a simple structure, e.g., which have a few arbitrary elements of one or two 
same arguments or whose equivalence groups are finite-dimensional. A number of results on 
group classification problems investigated within the framework of this approach are collected 
in [TUHUES] an d other books on the subject. 

Recently, based on the Lie-Ovsyannikov method and the investigation of the specific compat- 
ibility of classifying conditions, Nikitin and Popovych [65] developed an effective tool (we refer it 
as method of furcate split) for solving the group classification problem of nonlinear Schrodinger 
equation. In 2004, Popovych and Ivanova extended the method to complete group classifica- 
tion of nonlinear diffusion-convection equations by further considering the so called additional 
and conditional equivalence transformations [75]. In 2007, Ivanova, Popovych and Sophocleous 
present the extended and generalized equivalence transformation group, gauging of arbitrary el- 
ements by equivalence transformations for further investigation of nonlinear diffusion-convection 
equations [47]. Furthermore, Popovych and Ivanova et.al also extended these new group classi- 
fication idea to the nonclassical symmetries [571 190] and conservation law classification [49^179] . 
Up to now, these methods and different notations have been applied to investigating a number 
of different symmetry group, nonclassical symmetry and conservation law classification prob- 
lems paiaaiagiEiiisaiBaEaEMiiHiiisaisiiES!. 

However, almost all the research was concentrated on parabolic type nonlinear diffusion- 
convection equations and few of hyperbolic type nonlinear partial differential [37-39]. There- 
fore, the present paper is one of new extension of the above mentioned method and different 
notations to this classes of equations. The results of symmetry group, nonclassical symmetry 
and conservation law classification of class (JT[) present in this work are new. Hence, these will 
lead to some explicit applications in Physics and Engineering. 

The structure of the paper is as follows: 

Some physical examples contained in class (pQ) is discussed in section [2] 
In section [3l the complete group of usual equivalence transformations for class ([I]) and the 
extended one including transformations which are nonlocal with respect to arbitrary elements 
are constructed by using Lie-Ovsiannikov method and direct method. Taking into account the 
non-trivial subgroup of gauge equivalence transformations, we strongly simplify the solving the 
group classification problem to equation ([1]) with two different gauges g = 1 and g = h. 

The results of group classification under the extended equivalence transformation group for 
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the gauge g = 1 are contained in section 14,11 Then in section 14. 2\ the classification of gauge 
g = h are presented. The sketch of the proof of the obtained results are given in Section 14.31 
Classification with respect to the set of point transformations are presented in 14.41 We note 
that for both gauges two essentially different classifications are presented: the classification with 
respect to the (extended) equivalence group and the classification with respect to all possible 
point transformations. 

In Section exact solutions of some classification models are given by using the method of 
classical Lie reduction. 

After making a brief review of notation of singular reduction operator, we then carry out a 
preliminary analysis of nonclassical symmetry of the class ([T]) with the gauge g = 1 in section 
[6J As an example, we also present several reduction operators of a special nonlinear telegraph 
equation and constructed some non-Lie exact solutions for them. 

In Section the local conservation laws of these equations are exhaustively described. Using 
the most direct method, two classifications of local conservation laws up to equivalence relations 
are generated by both usual and extended equivalence groups. Equivalence with respect to these 
groups and correct choice of gauge coefficients of equations play the major role for simple and 
clear formulation of the final results. 

Finally, some conclusion and discussion are given in Section 

In the Appendix, classification results for the gauge g = 1 under the usual equivalence 
transformation group can be found. 

2 Physical examples 

Class (pQ) is a unified form of many significant second-order hyperbolic type nonlinear PDEs 
in Physics, Mechanics and Engineering Science. Physical examples corresponding to the case 
f(x) = g(x) = h(x) = 1 and K(u) = are collected in the well known paper [3|6], which 
describe the flow of one-dimensional gas, longitudinal wave propagation on a moving threadline 
and dynamics of a finite nonlinear string and so on. In what follows, we review several important 
physical models related with the coefficient f(x) ^ or h{x)K{u) ^ [88j . 

Example 1. Two-conductor transmission lines telegraph equation. The waves in two- 
conductor transmission lines having small transverse dimensions (in comparison with the char- 
acteristic wavelength) can often be described by the telegraph equation [53] 




where t is a spatial variable and x is time; I ,U, F(U), and G(U) are respectively the current 
in the conductors, the voltage between the conductors, the leakage current per unit length, and 
the differential capacitance. The form of F(U), and G(U) depend both on the configuration of 
the conductors, and on the properties of the medium filling it. 
Setting U = u, we obtain the telegraph equation 



which fall into ([TJ). 

Example 2. Longitudinal vibrations of elastic and non-homogeneous taut strings or bars. 
Suppose a string is taut along the x— axis. The equation giving the balance of momentum is 




I X = U U I t = F(U)U X + G(U) 



(4) 




(5) 



putt 



(6) 



where x is the coordinate of the point P in the present reference system, and 



x = x(y,t), 
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where y represents the coordinate of the corresponding point Pq of P in the reference shape, 
where uj = x — y, p is mass per unit length and T is the tension. To equation ([6]) we associate 
the following constitutive relations already considered in 



T = T(uj x ),p = p(x). (7) 

The balance law ([6]), with ([7]), transforms to the following second order partial differential 
equation 

T'( u ) . 

utt = [ -, -s u x \ x , . (8) 
pyx) 

where uj x = u. Of course ([8]) is particular case of (pQ). 

Example 3. Bar with variable cross section. The equation of motion of a hyperelastic 
homogeneous bar, whose cross sectional area is variable along the bar, is |22j 



where p is the (constant) mass density, uj = y — x is the displacement, y is the coordinate of the 
point P in the present reference system, x represents the coordinate of the corresponding point 
Pq of P in the reference frame, T is the tension and S{x) is the cross sectional area. 

Taking into account the constitutive relation T = T{uj x ) |22j . the equation @ becomes 

r T'(u) , ,S'(x). T(u) S'(x)T'(u) , , 

where u = u x . Obviously when ( )x = the equation (fTQ|) is included in JT]). 

Example 4. One dimensional propagation of visco-elastic stress waves. In |91j . the following 
constitutive laws were adopted for a nonlinear homogeneous visco-elastic model 

e t = $(T)T t + n(T), (11) 

with 

where T is the stress, uj is the displacement, do > 0, k and cq are constants. When the one- 
dimensional propagation of visco-elastic stress waves is investigated, combining the momentum 
equation 

putt = T x (13) 
with (fTT]) and taking into account (fT2j) . we obtain 

Txx = \l + d kc Q Tf Tt]t + (l + ZvT)^- (14) 
So, we fall into ([T]) when x and t are exchanged. 

Example 5. Hyperbolic heat equation. Many models for the heat propagation with finite 
speed give a hyperbolic equation which can be reduced to ([T]). In fact, quite recently, in order 
to describe one-dimensional heat conduction, the following partial differential equation has been 
considered in 184 



q Q a 2 5 2 2e 2 I 2q a 2 6 2 2 _ 

— + e > °™ ~ JT79 6t + r t ^ i e + ' x ~ ' (15) 
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where 9 denote the (absolute) temperature while qQ,a,^*,5 and e are suitable constants. The 
equation (|15p is based on a nonlinear model with relaxation. 
If we set 9 = 1/u, we obtain 

7* 1 
' r«t]H o~7t — ; — y>u t . (16) 



qoa 2 (5u + e) 2 qoa 2 (5u + e) 

This equation was obtained in [21] in the framework of Miiller's theory for heat propagation 
in rigid bodies. 

Example 6. One dimensional heat propagation in a rigid body. The models like Cattaneo's 
|18j . describing one dimensional heat propagation in a rigid body, are governed by a nonlinear 
equation of the type [25] 

o xx = [-c(9)e t ] t + -[[ c{e)de\ e e t . (17) 

where C{9) is the special heat, tq is the thermal relaxation time and x is the thermal conductivity. 
The equation (|16p and (|17p fall into (H|) when x and t are exchanged. 



3 Equivalence transformations and choice of investigated class 

In order to perform group classification of class (PQ) , we should first find its group of equivalence 
transformations. The usual equivalence group G~ of class ([1]) is formed by the nondegenerate 
point transformations in the space of (i, x, u, /, g, h, H, K), which are projectible on the space of 
(t,x,u), i.e. they have the form 

(t,x,u) = (T t ,T*,T u )(t,x,u), 

(f,g,h,H,K) = (Tf,T9,T h ,T H ,T K )(t,x,u,f,g,h,H,K), 

and transform any equation from class ([I]) for the function u = u(t,x) with the arbitrary 
elements (/, g, h, H, K) to an equation from the same class for the function u = u(t, x) with the 
new arbitrary elements (f,g,h,H,K). To find the connected component of the unity of 
we have to investigate Lie symmetries of the system that consists of equation ([1]) and some 
additional conditions, that is to say we must seek for an operator of the G~ in the form 

X = rd t + £d x + rjd u + ird f + ^d g + <t>d h + pd H + 9d K (18) 

from the invariance criterion of the following system: 

f(x)u tt = (g(x)H(u)u x ) x + h(x)K(u)u x , 

ft = fu = 0, 9t = 9u = 0, h t = h u = 0, H t = H x = 0, K t = K x = 0. 

Here u, f, g, h, H and K are considered as differential variables: u on the space (t, x) and /, g, 
h, H, K on the extended space (t,x,u). The coordinates r, £, r\ of the operator (fT8|) are sought 
as functions of t, x, u while the coordinates tt, <p, (ft, p and 9 are sought as functions of t, x, u, 
/, g, h, H, K. 

The invariance criterion of system (fT9l) yields the following determining equations for r, £, 77, 
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7r, p>, (j), p and 8: 

t x = r u = £ t = 6« = Vx = Vt = 0, T tt = Vuu = 0, 

TTt = VT« = VT/f = K K , (ft = ip u = ip H = <PK = 0, 

4>t = 4>u = 4>H = 4>K = 0, pt = Px = Pu = Pf = Pg = Ph = PK = 0, 

9 t = 9 x = 9f = e g = e h = 0, (20) 

7T CP 7T ip p 

+ 2£ x -2 n -^-pH = 0, - + 2Cx-2T t -^- = j-, 
f 9 f 9 H 

7T 

9xP + [-gixx + (2Tt - 6k — j )ftc + C^x + /zV 3 / + feCOg + fyz^ft - £x9x)H 

+h9 + (2r t - & - j + j)hK = 0. 

After easy calculations from (I20p . we can find the connected component of the unity of G~ for 
the class (HJ. 

Theorem 1. The usual equivalence transformation group G~ for the class ([I]) consists of the 
transformations 

t = 5\t + 5 2 , x = X(x), u = 8311 + 64, 

f = ^~/> 9 = e i e 2 lx x9, h = eieg l h, H = e 2 H, K = e 3 K, 

where 5j (j = 1,...,4) and 6j (i = 1,2,3) are arbitrary constants, 5i £3616263 ^ 0, X is an 
arbitrary smooth function of x, X x 7^ 0. 

It is shown that class (pQ) admits other equivalence transformations which do not belong to G~ 
and form, together with usual equivalence transformations, an extended equivalence group. We 
demand for these transformations to be point with respect to (t,x,u). The explicit form of the 
new arbitrary elements {f ,g,h,H,K) is determined via (t,x,u, f,g,h,H,K) in some nonfixed 
(possibly, nonlocal) way. We can construct the complete (in this sense) extended equivalence 
group G~ of class (PQ) by using the direct method. 

Theorem 2. The extended equivalence transformation group Gr for the class ([!]) is formed by 
the transformations 

t = 5\t + 5 2 , x = X(x), u = 6311 + 84, 

f = ^Y^f f e~ e4 f~i>dx, g = eie^Xxg f e~ e4 ^dx, h = €\e^ l h f e~ C4 ^dx, 
H = e 2 H, K = e 3 (K + e 4 H), 



where 5j (j = 1,...,4) and ei (i = 1,...,4) are arbitrary constants, 5i^3eie2e3 ^ 0, X is an 
arbitrary smooth function of x, X x 7^ 0, f~=J Jffidx. 

Remark 1. It should be noted that the existence of such equivalence transformations can be 
explained in many respects by features of representation of equations in the form ([I]). This 
form usually leads to an ambiguity because the same equation has an infinite series of different 
representations. In fact, two representations ^Q) with the arbitrary element tuples (/, g, h, H, K) 
and (/, g, h, H, K) determine the same equation if and only if 

f = ^rf f e ~ £i ^ dx > 9 = ei£2 1 X x g f e~ e4 ^dx, h = e^^h f e~ 64 ^ dx, 
H = e 2 H, K = e 3 (K + e 4 H), 

where 5i and ei (i = 1, . . . , 4) are arbitrary constants, €16263 7^ (the variables t, x and u do not 
transform!). The transformations (|2ip act only on arbitrary elements and do not really change 
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equations. In general, transformations of such type can be considered as trivial |63j ("gauge") 
equivalence transformations and form the "gauge" (normal) subgroup G~ 9 of the extended 
equivalence group Cr . Application of "gauge" equivalence transformations is equivalent to 
rewrite equations in another form. In spite of really equivalence transformations, their role in 
group classification comes not to choice of representatives in equivalence classes but to choice of 
form of these representatives. 

Remark 2. We note that transformations (|2ip with 64 7^ are nonlocal with respect to arbitrary 
elements, otherwise they belong to G~ and form the "gauge" (normal) subgroup G~ 9 of the 
equivalence group G~ . The factor-group G~/G~ 9 for class ([I]) coincides with G~/G~ 9 and can 
be assumed to consist of the transformations 

t = 5it + 52, x = X(x), u = 5 3 u + 5 4: , 

f~=if> ~9 = X x g, h = h,H = H, K = K, 

where Sj (j = 1, ... ,4) are arbitrary constants, 5\5 3 7^ 0, X is an arbitrary smooth function of 
x, X x + 0. 

Based on the this idea, we can gauge the parameter- function g in equation JT]) to 1. More 
exactly, using theorem [H we deduce that the transformation 

dx 

t = t,x= I , r ,u = u 



g{x) ' 

from G~ /G~ 9 reduce equation (P) to 

f{x)u- t i = (H(u)u x ) x + h(x)K(u)u s , 

where f(x) = g(x)f(x),g(x) = 1 and h(x) = h(x). (Likewise any equation of form (TjQ) can be 
reduced to the same form with f(x) = 1). That is why, without loss of generality we can restrict 
ourselves to investigation of the equation 

f{x)u tt = (H(u)u x ) x + h(x)K{u)u x , (23) 

Below, we denote and G± as the usual and extended equivalence transformation group of 
equation ([23]) . 



Theorem 3. The usual equivalence transformation group Gi for class (123jl consists of the 
transformations 

i = e 4 t + ei, x = e 5 x + e 2 , n = e 6 n + e 3 , 

f = ejes 2 e 7 f, h = e% l h, H = e 7 H, K = e^e 7 e 8 K, 

where e\, . . . e 8 are arbitrary constants and e^e^eee-jes 7^ 0. 



Note that for class (J23J) there also exists a non-trivial group of discrete equivalence transfor- 
mations generated by four involutive transformations of alternating sign in the sets {£}, {x, K}, 
{u}, {f,H,K} and {h,K}. Class (|23p admits other equivalence transformations being nonlocal 
with respect to arbitrary elements which do not belong to G^. We demand for these transfor- 
mations to be point with respect to (t,x,u). In such way, using the direct method we can find 
a generalized equivalence group Cr of class (P) . 

Theorem 4. The extended equivalence transformation group G^ for class (j23|) is formed by the 
transformations 



t = e 4 t + ei, x = e 5 / e eg f h dx + e 2 , u = e 6 u + € 3 , 

f = e 2 4 e^e 7 f fe- 2 ^f h dx, h = e 7 e^h J e"* S h da, H = e 5 e 7 H, K = e 8 (K + e 9 H), 
where e±, . . . eg are arbitrary constants and £163656769 7^ 0, Jh = Jh(x) dx. 



(24) 



9 



The group £r~ is a subgroup of Gr . Its transformation can be considered as from Gr , 
which preserves the condition g = 1. The transformations (|24p with non- vanishing values of 
the parameter eg are also nonlocal in the arbitrary element h. There exists a way to avoid 
operations with nonlocal equivalence transformations. More exactly, we can assume that the 
parameter-function K is determined up to an additive term proportional to H and subtract such 
term from K before applying equivalence transformations (122D . 

At the same time, there is another possible generalization of the gauge g = 1 to the general 
case of h, namely the gauge g = h. Any equation of the form ([T]) can be reduced to the equation 

f(x)u tt = {h(x)H(u)u x ) x + h{x)K(u)u x . (25) 
by the transformation t = t, x = j jnQdx, u = u from G~ /G~ 9 . 

Theorem 5. The usual equivalence transformation group for class (125j) consists of the 
transformations 

t = eit + e 2 , x = e 5 x + e 6 , u = e^u + ei, 

f = ej e^ 2 e 7 e 8 /, h = e s h, H = e 7 H, K = e^e 7 K, 

where e\, . . . e$ are arbitrary constants and e^ceeies ^ 0. 

Theorem 6. The extended equivalence transformation group G^ for class f)25f) is formed by the 
transformations 

t = e 1 t + e 2 , x = e 5 x + e 6 , u = e 3 u + e 4 , 

/ = efe 5 - 1 e 9 /e^ ; h = e 9 e^he^, H = e 5 H, K = e 7 (K - e 8 H), [ } 

where t\, . . . eg are arbitrary constants and e\e^e^e 7 e 9 / 0, fh = fh(x) dx. 
Remark 3. If H = 0, we assume h = 1 for determinacy. 

4 Group Classification of nonlinear telegraph equations 

In this section, we will present the group classification for the class (pQ) with the gauges g = 1 
and g = h under the extended equivalence transformations group Gr . 

Following the algorithm in [3,4,47,73,78] we are looking for an infinitesimal operator in the 
form 

Q = r(t, x, u)d t + x, u)d x + rj(t, x, u)d u (27) 

which corresponds to a one-parameter Lie group of local transformation and keep the equa- 
tion ([TJ invariant. The classical infinitesimal Lie invariance criterion for equation ([T]) to be 
invariant with respect to the operator (|27p read as 

pr^ Q(A) | A =o= 0, A = f{x)u tt - {g{x)H{u)u x ) x - h{x)K{u)u x . (28) 

Here pr^ 2 ) Q is the usual second order prolongation (TQlEi] of the operator (|271) . Substituting 
the coefficients of pr^ 2 ) Q into (|28[) yields the following determining equations for r, £ and r\: 

T~x — Tu — — — Vuu — 0, 

H(gn x ) x + hKrj x - fr] tt = 0, 

if Hi - fH£ - 2r t H - nH u + 2H£ X = 0, 

i.9xV + 2gr] x )H u + [(2r] xu - £ xx )g + (2r t -t[ x - £j-)g x + i9xx]H + hrjK u (29) 
+(2r t - Cx - ^ + Zfy)hK = 0, 



2rj tu -T tt = 

2(6. - Vu)H u - 2r t H u - nH uu + ^H u £ - fH u £ + 7] U H U = 0. 
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Investigating the compatibility of system (|29p we find that the final equation of system (|29p 
is an identity (substituting the third equation of system (|29[) to the final one can yield this 
conclusion). With this condition, system (|29p can be rewritten in the form 

r x = T u = it = iu = Vuu = 0, 2i] tu = T tt , (30) 

2fe - Ti) + ( f -f - = ^fV, (31) 

(g Vx ) x H + - frjtt = 0, (32) 

fx 

(g x rj + 2gr] x )H u + [(2^ - + (2r t - - £y )ftr + £g X x]H + frr/i^ 
+ (2T t -&-£^+£-^)/iiT = 0. ( 33 ) 
Equations (|3Up do not contain arbitrary elements. Integration of them yields 

r = r(t), Z = £(x), r, = vHt,x)u + V°(t,x), r]\t, x) = -r t + a(x). (34) 

Thus, group classification of (JH) reduces to solving classifying conditions (|3T1) ~ ([33]) . 

Splitting system (|3ip --(|33 p with respect to the arbitrary elements and their non-vanishing 
derivatives gives the equations = 0, ^ = 0, i] = for the coefficients of the operators from yl ker 
of (fT]). As a result, we obtain the following assertion. 



4.1 Classification under the gauge g = 1 

Theorem 7. The Lie algebra of the kernel of principal groups of ([1]) with the gauge g = 1 is 

A ker = (d t ). 

Theorem 8. A complete set of inequivalent equations ([I]) with the gauge g = 1 with respect to 
the transformations from G± with A ma,x ^ A kcr is exhausted by cases given in tables[I^\3i 



Table 1. Case of \/H(u) (gauge g — 1) 



N 


K(u) 




h(x) 


Basis of A max 


1 


V 


V 


V 


d t 


2a 


V 


e px 


1 


d t , ptdt + 2d x 


2a' 


1 


\x\" 


x- 1 


dt, (p + 2)tdt + 2xd x 


3 


1 


x- 2 


x- 1 


dt,td t + xd x 


4 





1 


1 


dt, d x ,tdt + xdx 



Here p £ {0, 1} mod G~ in case [21 p 7^ —2 caselJK'. 
Additional equivalence transformations: 

1- Hi (p = 0, K = 1) ->EK (p = 0, K = 0): t = t,x = x + t,u = u; 
2. Hi' (p = -2, if = 1) -> UK' (p = -2, JC = 0): t = t, 5 = a;e*, w = u. 



Table 2. Case of H(it) = e" u (gauge g = 1) 



N 


K{u) 


f(x) 


h(x) 


Basis of A max 


1 





V 


1 


d t , tdt — 2d u 


2 




\x\ p 


\x\ q 


d t , \p(u -v) + (u- 2v) - qfi]td t + 2(ji - v)d x + 2(q + l)d u 


2* 




e px 


ee qx 


d t , ]p(ji -v)- q^]td t + 2(u - v)d s + 2qd u 


3 


ue u 


h 2 e"f h 


(h- 1 )" = -2ph 


dt, (2p + q)td t + 2h- l d x - 4pd u 


4 


e uu 


1 


1 


d t ,d x , (u — 2v)td t + 2{u — v)xd x + 2d u 


5 







1 


dt,td t - 2d u ,atd t + 2{/3x 2 + j lX + jo)d x + 2/3xd u 


6a 





1 


1 


d t ,d x , tdt - 2d u , tdt + xd x 


6b 





x- 3 


1 


d t , xd x — d u ,x 2 d x + xd u , tdt ~ 2d u 
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Here (p, v) £ {(0, 1), (1, v)}, v 7^ p in cases [5J [2f and 2] /1 = 1 and v 7^ 1 in the other cases; g 7^ — 1 in case [2T 
(otherwise it is subcase of the case 1.2a); e = ±1 in case^; a, ,5,71,70 =const and 

fHx) = exp / [ ^l—^l±^dx 

Additional equivalence transformations: 
1. [6]b — > [6^1 : t = tsignx, x — 1/x, u = u — ln\x\. 



Table 3. Case of H(u) = w M (gauge g = 1) 



N 




#(«) 


f(x) 


ft(a;) 


Basis of A max 


1 







V 




1 


dt, ptdt — 2ud u 


2 


V 




|af 






d t , [{p - q)p — pu + fi — 2v]td t 














+2(p - v)xd x + 2(q + l)ud u 


2* 


V 


i«r 


e px 






dt, [(p-q)p-qv]d t 














+2(p - v)d x + 2qud u 


3 


V 


liiC In |m| 


/j 2 e 9 / hdx 


(ft" 1 


)" = -2pft 


dt, (2pp + q)td t + 2ft" 1 dx - 4pud u 


4 





V 


h 2 


(ft- 


-i)" = 


dt, h~ 1 d x 


5 





u 


ftV h 


(ft- 1 


)" = -2pft 


dt,td t + 2h- 1 dx-4pd u 


6 


V 


M" 


1 




1 


d t , d x , (p - 2u)td t + 2{n - u)xd x + 2ud u 


7 


^-4 





/ 3 (z) 




1 


d t ,ptd t — 2ud u , 














atdt + 2[(p + l)f3x 2 + jix + jo]dx + 2/3xud u 


8a 


^ -4 





1 




1 


d t ,ptd t — 2ud u , d x ,tdt + xd x 


8b 


7^ -4 -± -1 









1 


dt,ptd t - 2udu, {n + 2)tdt - 2{p + l)xdx, 














(p + l)x 2 d t + xud u 


8c 


-1 





e* 




1 


dt , tdt + 2udu , dx — ud u , 














td t + xd x — xudu 


9 


-4 





/ 3 (z)U=-4 




1 


dt, 2tdt + ud u ,t 2 dt + tudu 














atdt + 2{-3/3x 2 + 7 ix + 70)^ + 2/3xud u 


10 


-4 





1 




1 


dt, 2tdt + ud u ,dx,t 2 dt + tudu, 














2xdx — udu 


11 


4 
3 





1 




1 


dt , 2tdt + 3ud u , dx , tdt + xdx , 














x 2 d x — 3xud u 



Here v 7^ p; e — ±1 in case[2] (otherwise it is subcase of the case 1.2a); a, /3, 71, 70 =const, and 



/ (x) = exp 



-(3At + 4)ftx-27i +q 
(p + l)/3x 2 + j±x + 70 



dx 



Additional equivalence transformations: 

1. Eb-^EJi: t = t,x = -l/x,u = \x\~ r ^u; 

2. [8t -> [HJi (p = -1): t = t, x = x, u = e x u. 

Remark 4. Tables [T]{3] are the results of classification for class ([1]) with the gauge g = 1 with 
respect to the extended equivalence transformations group GJ". The classification for class (JTJ) 
with the gauge g = 1 with respect to the usual equivalence transformations group Gi are very 
complicated. In some cases, the determining equations can not be solved explicitly. Therefore, 
We list these results as an Appendix. 

Remark 5. The proof of theorem [8] follows directly from the analysis of the Section f4. 31 
4.2 Classification under the gauge g = h 

Theorem 9. The Lie algebra of the kernel of principal groups of ([1]) with the gauge g = h is 

A^ = {d t ). 
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Theorem 10. A complete set of inequivalent equations ([T]) with the gauge g = h with respect 
to the transformations from G7~ with A max ^ A ker is exhausted by cases given in tables [JEJ 



Table 4. Case of VH(u) (gauge g — ft) 



N 


K{u) 




h(x) 


Basis of A max 


1 


V 


V 


V 


dt 


2a 


V 


e px 


1 


d t , ptd t + 2d x 


3 


1 


x- 2 


x- 1 


d t ,tdt + xd x 


4 





1 


1 


d t , d x ,td t + xd x 



Here p G {0, 1} mod in case [2] 



Table 5. Case of H[u) = e" u (gauge g = h) 



N 


K{u) 


/(*) 


h{x) 


Basis of A max 


1 





V 


1 


d t , td t — 2d u 


2 


e" u 


\x\ p 


\x\« 


d t , [(p-q+l)fj,-(p-q + 2)u]td t + 2{p - v)xd x + 2d u 


3 


ue u 


e px 2 + qx 


e px 2 


d t , {2p + q)td t + 2d x ~ 4pd u 


4 


e vu 


1 


1 


d t ,d x , (/J, — 2v)td t + 2(fi — v)xdx + 2d u 


5 







1 


d t ,tdt - 2d u , octdt + 2(/3x 2 + ^x + j Q )d x + 2/3xd u 


6a 





1 


1 


dt,d x , tdt — 2d u , tdt + xd x 


6b 





x- 3 


1 


d t , xdx — d u , x 2 d x + xd u , tdt — 2d u 



Here (fj., v) G {(0, 1), (1, v)}, v 7^ ji in cases [2] and [4] [i — 1 and v 7^ 1 in the other cases. 



Table 6. Case of H(u) = tt M (gauge g = ft) 



N 






m 




Basis of A max 


1 







V 


1 


dt, utd t — 2ud u 


2 


V 


\u\ u 




m« 


d t , [{p-q + l)ii-(p-q + 2)v\tdt 












+2(p - v)xdx + 2ud u 


3 


V 


\u M ln m 


e px 2 +qx 




d t , (2p(i + q)td t + 2d x — 4pud u 


4 





it 


ftV h 


(ft- 1 )" = -2ph 


dt,tdt + 2h- 1 d x -'ipd u 


5 





V 


h 2 


(ft" 1 )" = 


d t , h~ 1 d x 


6 


V 


M M 


1 


1 


d t , d x ,(u — 2v)td t + 2(fi — v)xdx + 2ud u 


7 


^-4 





/ 3 (z) 


1 


dt,fJ,td t — 2ud u , 












atdt + 2[(p + l)/3x 2 + 71a; + -y ]d x + 2f3xud u 


8a 


^ -4 -4 





1 


1 


dt,/J,td t — 2ud u , d x ,tdt + xd x 


8b 


7^ -4 -1 
/- ^> 35 1 







1 


d t ,utdt — 2ud u , (u + 2)td t — 2(/i + l)xd x , 












(fi + l)x 2 d t + xud u 


8c 


-1 







1 


dt,td t + 2ud u , d x — ud u , 












tdt + xdx — xudu 


9 


-4 





/ 3 («)U=-4 


1 


dt, 2td t + ud u ,t 2 d t + tud u 












atdt + 2{-3/3x 2 + 7 ix + 70)^ + 2/3xud u 


10 


-4 





1 


1 


d t , 2td t + ud u ,dx,t 2 d t + tud u , 












2xdx ~ ud u 


11 


4 
3 





1 


1 


dt , 2td t + 3ud u , d x , tdt + xdx , 












x 2 dx — Sxudu 



Here v =f= /*. 



In tables HH3] and HH6] we list all possible G~ -inequivalent sets of functions f(x), h(x), H(u), 
K{u) and corresponding invariance algebras under the gauges g = 1 and g = h respectively. We 
give the same numbers for the corresponding (G~-equivalent) cases in the gauges g = 1 and 
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g = h. The asterisked cases from tables [2] and [3] are equivalent to the cases from tables [5] and 
[6] with the same numbers, where the parameter-function h takes the value h = x. The similar 
non-asterisked cases correspond to the same cases from tables [5] and [H where p' = jSfj?' = — +rr 

orp' = -£p2 if g = p + l. 

In what follows, for convenience we use double numeration T.N of classification cases, where 
T denotes the number of the table and iV the number of the case in table T. The notation 
'equation T.N' is used for the equation of the form ([!]) where the parameter functions take the 
values from the corresponding case. 

The operators from tables [TH3] or HHS form bases of the maximal invariance algebras if the 
corresponding sets of the functions /, h, H, K are G~-inequivalent to ones with most extensive 
invariance algebras. For example, in case 13111 the adduced operators have the above property iff 
f + f- 

Remark 6. Case 1112a is equivalent to case 1112a ' with respect to transformation t = t,x = 
ln\x\,u = u,H = H,K = K — H,p = p + 2 from G~. We adduce case 1112k ' here for the 
convenience of presentation of results only. 

4.3 Proof of classification results 

Now, let us use the method of furcate split [47)65,78] to prove the main classification theorems [8] 
and [lOl It should be noted that it seems impossible to formulate complete results of group 
classification of class ([I]) with respect to usual equivalence group G~ in a closed form. This can 
be seen from the classifications for the gauge g = 1 adduced in the Appendix, while it is quite 
easy to solve the problem of group classification with respect to the extended equivalence group 

cr. 

The basic idea of the method is based on the fact that the substitution of the coefficients 
of any operator from the extension of A ker into the classifying equations results in nonidentity 
equations for arbitrary elements (see [UlEUEE] ^ or more details about the method). In the 
problem under consideration, the procedure of looking for the possible cases mostly depends 
on equation pip . For any operator Q £ ^4 max equation (I3ip gives some equations on H of the 
general form 

(au + b)H u = cH, 

where a, b, c are constant. In general, for all operators from A max the number k of such 
independent equations is not greater than 2; otherwise they form an incompatible system on H. 
k is an invariant value for the transformations from Cr . Therefore, there exist three inequivalent 
cases for the value of k: 

1. k = : H(u) is arbitrary, 

2. k = 1 : H(u) = e^ u or H(u) = Qji / 0) mod GT, 

3. k = 2: H{u) = 1 mod G~. 

Furthermore, in order to provide the final presentation of classification results in a simple way, 
the choice of a gauge for the arbitrary elements is very important for solving the determining 
equations. It is more convenient to constrain the parameter- function g instead of / in class (pQ). 
The next problem is the choice between gauges of g. The case Ke(1,H) and k > 1 is easier to 
be investigated in the gauge g = h. In the other cases we obtain results in a simpler explicit 
form and in an easier way using the gauge g = 1. Let us consider these possibilities in more 
detail, omitting cumbersome calculations. 
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Case 1: k = (the gauges g = 1 and g = h, tables [T] and EJ). We first consider the case 
gauge 5 = 1. Since H{u) arbitrary, this means that the coefficients of any operator from A max 
must satisfy r\ = and 

2fe - r t ) + = 0, (35) 
-if (£h) x + = 0. (36) 

(i) Let us suppose that K g" (1,H). It follows from equation (i36j) that £ x = 0. Therefore, 
equation ([35]) must be in the form / x = without fail. Solving this equation yields cases [2^,. 

(ii) Now let K G (1,H), i.e. if = 5 mod where 5 G {0,1}. Then equation (j36j) can be 
decomposed into the following ones 

£** = 0, 8(£h) x = 0. (37) 

Integrating of the latter equations up to G^ results to cases (2^'-[3] of table [TJ 
The classifications for the gauge g = h can be derived in a similar way. 

Case 2: k = 1 (the gauges g = 1 and 5 = h, tables El [3] and U M ) ■ Here H G {ef M ,vt t ,fi^ 0} 
mod Cr and there exists Q £ J 4 max with 77 7^ 0, otherwise there is no additional extension of 
the maximal Lie invariance algebra in comparison with the case k = 0. Below we consider the 
gauge g = h in details, the gauge g = 1 can be proved in a similar way. If H = e^ n we assume 
H = l. 

Case 2.1: Let us investigate the first possibility H = e M " (table [5]). Equations ([3T]) and ([34"|) 
imply r/ w = 0, i.e. 77 = rj l (t,x) and r tt = 0. Therefore, equation ([33]) looks like K u = vK + Aii 
with respect to K, where is, b = const, otherwise rj = 0. 

Consider first the case Kg(1,H). Under the above suppositions, equations (i3~l~T) - (j33~j) can be 
rewritten as 

^t = {2v-n)r ] 1 + 2T U (38) 
vl = ril = 0, i xx = T tt = 0, (39) 
^ = (M-^y, (^), = -Ar? 1 . (40) 

Here and below (p = f /h. From equation (f38|) we can get 99 G {e 9X , |x| r (r 7^ 0), 1} mod G%. 

For ip = e qx it follows from the determining equations (j38l) - (j40"j) that £ x = 0,v = fi, \ 
0, 2j» = 2a. Thus, /i = /i e a:r2+/li:!; = e Qx2 mod G%, a ^ 0J = hp = e ax2+qx ,K = Xue u mod 
G^ that falls precisely into case 15131 

Hp = \x\ r ,r / 0, thenr£/x = (2u-/j,)r] 1 +2r t . Therefore, £ = (^-v)t] 1 x, (fi-u)('^) x = -A. 
Since /i 7^ f (otherwise, K G (1, H)) we have A = mod G£. Therefore, /i = mod G~. 
Then / = |x| p ,p 7^ g, and we obtain case 15121 

Value if = 1 results in 2r t = (// — 2u)r] 1 , £ = (fi — v)r] l x + £o- \iv = [i then A 7^ (otherwise, 
if G (1,-ff)), A = 1 mod G£ , (fy) x = 2a(r ? 1 = -2a). Therefore, h = h e ax * +hlX = e ax2 
mod G^ , if = ue u that follows to case 15131 If v 7^ fj,, then A = mod G£ . Therefore, 
h G {l^l 9 , 1, e px } G^ that yields subcases of 15131 15l2l and case 15141 correspondingly 

Now, we consider the case K G (1,H), H T^const. In contrast to the previous case, it is 
more convenient to consider this case using the gauge g = 1. In such case K = 0, 1 mod G^ . 
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Application of the above suppositions reduces the determining equations to the system 

fx 

2£ x + -ji = fir] 1 + 2r t , vL = 0, Krj x = (prj tt , 
(£* + — £ - n)K = 0, U = W X - 

Note that rj l = ^r t +a(x), thus we have rj 1 = ^r t +f3x+ao and £ = /x/3x 2 +7ix+7o. Substituting 
these values into the first determining equation we obtain 

fx 1 
{iij3x 2 + 71X + 7o)y = -3/i/3x + -(/x + 4)r t + ria - 271. 

This equation gives I linearly independent equations for / of form (a 2 x 2 +a 1 x+a°)^j- = /3ix+/3o- 
If I = 0, then ^ = 0, j3 = 0, + 4)rt = — /xao- Considering case I = 1, we get (a 2 , /3i) / (0, 0), 
(a , /3b) / (0,0), otherwise / > 1. At last, if I > 2, then / G {l,e px , |x| p ,p ^ 0} mod G^. 

Direct substitution of the above values to the determining equation for K = and obvious 
integration leads to the cases I5TT1 (case / = 0), 15151 (case / = 1) and 1516a . 1516b (case 1 = 2). 

Classification in case K = 1 is more cumbersome, and corresponding results can be reduced 
to cases I5l3l I5l2l and [5l4l 

Case 2.2: Consider the case H = u' 1 (table [6]). Equations ([3T]) and (fM)) imply rj = (^74 + 
a(x))u = r] 1 (t,x)u. Therefore, equation (f33l) with respect to K looks like uK u = vK + XH, 
where v,b = const, otherwise rj = 0. 

Let Kg(1,H). Using the above suppositions, we can rewrite equations (|3T]) - (f33l) as 

= (2v - fir? + 2r t , 

nl = vl = 0, &x = T tt = 0, (41) 

Z x = (li-v)r,\ (^) X = -X V \ 

From the first equation of system ()41|) we can get <p G {e qx , |x| r (r 7^ 0), 1} mod G^. 

For p = e qx it follows from the determining equations (|4ip that £3 = 0, v = fi, A 7^ 0, = 2a. 
Thus, /i = h e ax2+hlX = e ax2 mod G% , a ^ 0, f = hip = e^+^if = X\u\^\n\u\ mod G£ that 
falls precisely into case 16131 

If ip = \x\ r ,r / 0, thenr^/x = (2i/- /x)?? 1 +2r t . Therefore, £ = {ii-v)r] l x, (fi-v){^) x = -A. 
Since /j / 1/ (otherwise, if G (1,H)) we have A = mod Therefore, h = \x\ q mod G~. 
Then / = |x| p ,p 7^ q, K = \u\ u , and we obtain case 16121 

Value <p = 1 results in 2r t = (ri— 2u)rj 1 ,Z = (ri — v)r] l x + ^o- If v = A 4 then A 7^ (otherwise, 
K G (1,-ff)), A = 1 mod G£ , (^ = 2a(r? 1 = -2a). Therefore, /i = /i e Qx2+ ' llx = e Qx2 
mod , ^ = I tt I ^* In I I that follows to case 16131 If v 7^ fj,, then A = mod G^_ . Therefore, 
h G {l^l 9 , 1, e px } G^ that yields subcases of 16131 16l6l and [Tl2l correspondingly. 

Now, we turn to the case K G (1,H), H T^const. For convenience we also consider this case 
using the gauge g = 1. Hence K = 0, 1 mod G^ . Application of the above suppositions reduces 
the determining equations to the system 

2£r + j £ = /V + 2r t , r4 = 0, if r/i = pril , 
(£x + — Z - n)K = 0, ^ = 2(^ + 1)^. 

Solving this system and noting that r/ 1 = \r t + a(x), we can get r/ 1 = + /3x + ao and 
£ = (/x + l)/3x 2 + 71S + 7q. Substituting these values into the first determining equation and 
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differentiating it with respect to the variable t we obtain 

(42) 



|(/i + 4)r tt = 0. 

((/i + l)/3x 2 + 7 ix + 70)^ = -(3/U + 4)/3x + ±(/i + 4)r t + /xa - 2 7 i- 



The first equation of the above system implies that there exist two cases should be considered: 
Tu = if fj, 7^ —4 and r« 7= if /z = — 4. The second equation of system (|42p gives Z linearly 
independent equations for / of form (a 2 x 2 + o^x + a°)^j- = fi\x + /?o- If Z = 0, then £ = 
0,/3 = 0, ^(/x + 4)t 4 = —ficto- Considering case I = 1, we get (a 2 , /?i) 7^ (0,0), (a°,/?o) 7^ (0,0), 
otherwise Z > 1. At last, if I > 2, then / £ {1, e^, / 0} mod G7\ 

For the case = 0,// ^ —4, substituting the above values into the determining equation 
directly for if = and obvious integration leads to the cases Kill I (case 1 = 0), 16171 (case I = 1) 
and 1618a . 1618b . 1618b . 161 1 1 1 (case 1 = 2). The case tu 7^ 0, /j, = —4 with K = is corresponding to 
the results EM I6TT01 

The classification for K = 1 is corresponding to subcases 16121 10131 and 16161 



Case 3: k = 2 (the gauges g = 1 and g = h, tables O [3] and E] ) . The assumption of two inde- 
pendent equations of form (|3ip for H yields H = const, i.e. H = 1 mod G~. -K^ 7= (otherwise, 
equation (pQ) is linear). In what follows, we only use the gauge g = h. Equations (|3T]) - ([33l) can 
be written as 

2fe-r,) + (^-^)e = 0, (43) 

(hr] x ) x + hKrj x - frju = 0, (44) 

K u r, + i x K + U + 2rjl + {i) x h x = 0. (45) 

The latter equation looks similar to {au + b)K u = cK + d with respect to K, where a,b,c,d = 
const. Therefore, to within transformations from G~, K must take one of four values: 

K = u u , v^Q,l, K = lnu, K = e u , K = u. 

Classification for these values is carried out in the way similar to the above. The obtained 
extensions can be entered in either table [5] or table [6j The gauge g = 1 can be proved in a 
similar way. 

The problem of the group classification of equation (pQ) is exhaustively solved. 

4.4 Classification with respect to the set of point transformations 

Although we have performed the classification by using extended equivalence group, we can find 
in the classification results equations that some cases from tables [THH] or 0HH] are equivalent with 
respect to point transformations which obviously do not belong to G~. These transformations 
are called additional equivalence transformations and lead to simplification of further application 
of group classification results (see reference (37JII7JI78] for details). The simplest way to find 
such additional equivalences between previously classified equations is based on the fact that 
equivalent equations have equivalent maximal Lie invariance algebras. Explicit formulas for pairs 
of point-equivalent extension cases and the corresponding additional equivalence transformations 
are adduced after the tables. One can check that there exist no other point transformations 
between the equations from tables [TH3]or tables HHEJ Using this we can formulate the following 
theorem. 

Theorem 11. Up to point transformations, a complete list of extensions of the maximal Lie 
invariance group of equations from class ([T]) is exhausted by the cases from tables [7H3 or tables 
fj| numbered with Arabic numbers without Roman letters and subcases "a" of each multi-case. 
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As one can see, the above additional equivalence transformations have multifarious structure. 
This displays a complexity of a structure of the set of admissible transformations. Usually the 
problems of finding of all possible admissible transformations are very difficult to solve, see, 
e.g., [52,55,76,80]. We will try to discuss the structure of the set of admissible transformations 
of class (UJ) in a sequel paper. 

5 Lie reduction and similarity solutions 

In this section new Lie exact solutions for the equations from the initial class are constructed 
to just illustrate possible applications of the classification results obtained. We mainly perform 
group analysis of three classes of equations possessing nontrivial symmetry properties from 
the obtained classification lists by the reduction method and then apply to finding similarity 
solutions. For this purpose, we first construct the optimal sets of subalgebras for each kind of 
maximal Lie invariance algebras arising from group classification, then perform the reductions 
with respect to obtained subalgebras. The method of reduction with respect to subalgebras of 
Lie invariance algebras is well-known and quite algorithmic to use in most cases; we refer to the 
standard textbooks on the subject [701173] . 

We first consider the case [6] of Table [U i.e., the equation 

u tt = (u^u x ) x + u u u x , (46) 
which admits the three-dimensional Lie invariance algebra q generated by the operators 

Qi=d t , Q2 = d x , Q 3 = (fx - 2v)td t + 2((j, - v)xd x + 2ud u . 
These operators satisfy the commutations relations 

[Ql, Q 2 ] = 0, [Qi, Q 3 ] = (m- 2i/)Qi, [Qi, Q 3 ] = 2(/i - u)Q 2 . 

An optimal set of subalgebras of the algebra q can be easily constructed with application of 
the standard technique |70U73j . Another way is to take the set from [73], where optimal 
sets of subalgebras are listed for all three- and four-dimensional algebras. A complete list 
of inequivalent one-dimensional subalgebras of the algebra g is exhausted by the subalgebras 
{Qi}, {Q2), {Q3), {Q2 — Qi), {Q2 + Qi)- This list can be reduced if we additionally use the 
discrete symmetry (t, x, v) — > (t, —x, v), which maps (Q2 + Qi) to (Q2 — Qi), thereby reducing 
the number of inequivalent subalgebras to four. 

The optimal set of two-dimensional subalgebras is formed by the subalgebras (Q3, Qi) , 
{Q3, Q2), and (Qi, Q2). Lie reduction to algebraic equations with the latter two two-dimensional 
subalgebra leads only to the trivial zero solution. Below we list all the other subalgebras from 
the optimal set as well as the corresponding ansatze and reduced equations in Table [7J Solutions 
of some reduced equations are adduced. 



Table 7. Reduced ODEs and algebraic equation for equation (146 \ . 



N 


Subalgebra 


Ansatz 


Reduced ODE 


1 

2 
3 

4 
5 


(Qi) 
(Q3) 

{Q2-Q1} 
(Qs,Qi) 


u = ifi(uj), Lj = X 
u = <fi(uj), Lj = t 
u = t a if(oj), w = Jr, 
a = % , p = 2 ^~ v) 

fi—2u ' " f± — 2v 

U — fl{oj), LO = X + t 
1 

u — Cx 


{<fP<Pu)u + ^tfiu = 

<Akw = 

a(a - l)<p + (/3 2 - 2a/3 + 2/3)cj<p u + /3 2 u 2 <p uu 
— (<P^<Pu)lj — <p v tptj = 

{ip^tpuj)^ + ^ifu, — = 

(1 + v)C^ +1 + (/lj - v)C v+1 = 
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Two kinds of explicit solutions can be constructed for arbitrary values of [i and v: the x-free 
solution u = Co + c\t and the stationary solution u = {jfzr-) v ~»x»- v . We can also construct two 
implicit solutions from the first and the fourth ODEs in table [7J 

U = if{x), u = ip(x + t), 

where cp and ip sarisfy 

' -<P^ U + — ^~r x + f ~%dx = 0, — !— + -^~ v + -^—u> + I -^rduj = 0, 



ix-v v + 1 J ip u+1 ' ii-v r V V + l J lfj u+1 

and u> = x + t, C\,C 2 are arbitrary constants. Let \x and v be particular values, we can get 
some number of explicit exact solutions from the above implicit solutions. For example, from 
the latter implicit solution, we can get four triangular function exact solutions if \i = 1, v = 2: 

u = -\ - ^tanp^(x + t)], 
u = -i + ^cot[^(x + i)], 

u = -§ - & tan[^(x + t)\ ± i^/l + 3tan 2 [^(x + t)], 
u = ~\ + cot[^-(x + t)] ± i^/l + 3cot 2 [^(x + t)]; 

and a rational solutions if // = — l,i> = — 2: 

u = x + t. 

Lie reduction and exact solutions of 'truly' variable-coefficient nonlinear telegraph waves are 
most interesting. We consider two cases I2l2f and 16121 i.e. equations 

e P x u tt = (u'"u x ) x + eeV Ux , (47) 



\x\ p u tt = {\x\ q u^u x ) x + \x\ q v?u x . (48) 

For each from these cases we denote the basis symmetry operators adduced in Table [2] and [6] 
by 0i = (Qi = d t , Q 2 = [(p - q)fi - qv]d t + 2(/i - v)d x + 2qud u ) and g 2 = (Q x = d t , Q 2 = 
[(p — q + l)fJ- — (p — q + 2)v]tdt + 2(/i — v)xd x + 2ud u ), which are all non-commutative algebra. 
A complete list of inequivalent non-zero subalgebras of X or cj 2 is exhausted by the algebras 
(Ql), (Q 2 ) and (Qi, Q 2 ). 

Lie reduction of the equations f4"T|) and (j!5| to ordinary differential equations (ODEs) and 
an algebraic equation can be respectively made with the one-dimensional subalgebra (Qi), (Q 2 ) 
and the two-dimensional subalgebra (Qi, Q 2 ) which coincides with the whole algebra X or fj 2 - 
The associated ansatzes and reduced equations are listed in Table M and EJ 



Table 8. Reduced ODEs and algebraic equation for equation (I47|l . 



N 


Subalgebra 


Ansatz 


Reduced ODE 


1 

2 
3 


(Qi) 
(Qi,Q2> 


u = (tp(w)) M+ 1 , UJ = X 

u = exp(y(w)), oj = a; 
u=\t\ a tp(u), u> = x + p]n\t\ 

a _ 2 9 o _ 2 P("-M) 

(<3-p)M+Pf' ^ P[(j>— P"] 

U = Ce i^-" 


¥W + ee* w ^<p^TT _1 = if H^-l 
^ u + ee<" +1 >* ,+ «> << , = if fi=-l 

a(a - l)tp + (2a/3 - /3)tp u + p 2 ip uul 
— lq~ v (ip^ip^u — euj q ~ p ip u if UJ = 

q(ji + 1)C M+1 + e(p - v)C v+1 = 



Table 9. Reduced ODEs and algebraic equation for equation (I48|) . 
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N 


Subalgebra 


Ansatz 


Reduced ODE 


1 


<Qi) 


u = 


(y>(w))»H-i, w 


= a; 


tpfuu — V 3 ^ + qffui + ip w <P t>+1 =0 if n=fc—l 






u = 




= a; 


+w«e<" +1 >*V w =0 if M = -l 


2 


(Q-2) 


u = 


t\ a (p(to), L0 = 




a (a - + (/3 2 + 2a/3 - /3)a;^ + P 2 U) 2 ip uu 






a = 
= 


2 




-u- p (wVV)« - = o 






(p-<?)(a»-")+(a«- 


2k) ' 












3 


<Qi,Q 2 ) 




u = Cxv-" 




[(g - 1)( P -v)+n+ 1]C» +1 v)C v+1 = 



Reduction to algebraic equations gives the following solutions of the initial equations (|47p and 
respectively: 



. q(M + i; 

[I — V 



1i 



(g-l)(/x -»/) + (// + !) 
v — fj, 



M+l 

w+1 1 



Furthermore, some of the reduced ordinary differential equations in tables [5] and are the 
modification of the Emden- Fowler and the Lane-Emden equations |35tl75|. For example, the 
first equation corresponding to case 1 of tableware the standard Emden- Fowler equation, while 
the second one to case 1 of table [9] is the generalized Lane-Emden equation. Solutions of these 
equations are known for a number of parameter values (see e.g. j35|.IT5J ) . As a result, classes 
of exact solutions can be constructed for wave equations (|47p and (|48p for a wide set of the 
parameters \x and q. We omit these results in order to avoid a cumbersome enumeration. 



6 On nonclassical symmetries 

In 1969, Bluman and Cole introduced an essential generalization of Lie symmetry in the study 
symmetry reduction of the linear heat equation [12J. These generalized symmetries are often 
called nonclassical symmetries (called also conditional or Q-conditional symmetries) nevertheless 
it was not used in [12] . A precise and rigorous definition of this notion was suggested noticeably 
later [30, 96J (see also [58] for a recent discussion on definition of nonclassical symmetries). 
Since then there is an explosion of research activity in the area of investigation of nonclassical 
symmetries of PDEs arising from different fields of physics, biology and chemistry |20 y21l[36"ll6H 
I67TI69] . Some of these works concern with nonlinear wave equations. See, for example, |27]I44]. 

Generally speaking, there are two main features of nonclassical symmetries of differential 
equation difference from Lie symmetries. The first one is that they can yield solutions not 
obtainable from the classical Lie symmetries. The second feature is the deriving systems of 
determining equations for nonclassical symmetries which crucially depends on the interplay be- 
tween the operators and the equations under consideration, and thus are different from the 
Lie symmetries. Due to these facts, when studying the general form of nonclassical symmetry 
operators Q = r(t,x,u)dt + x, u)d x + r](t, x, u)d u ((r, £) ^ (0,0)) for the linear heat equa- 
tion ut = u xx , there are two essentially different cases of nonclassical symmetries should be 
considered: the regular case r / and the singular case r = 0. The factorization up to the 
equivalence of operators gives the two respective cases for the further investigation: 1) r = 1 
and 2) r = 0, £ = 1. In particular, for the singular case the system of determining equations 
for nonclassical symmetries consists of a single (l+2)-dimensional nonlinear evolution equation 
for the unknown function r\ and, therefore, is not overdetermined. The determining equation is 
reduced by a nonlocal transformation to the initial equation with an additional implicit inde- 
pendent variable which can be assumed as a parameter [29]. The linearity of the heat equation 
is inessential here. 

Recently, based on the above discoveries, Kunzinger and Popovych [57] raise a number of 
interesting questions, to wit: Is the partition of sets of nonclassical symmetry operators of the 
linear heat equation with the conditions of vanishing and nonvanishing coefficients of operators 
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to regular and singular cases universal for any differential equations or is it appropriate only 
for certain classes of differential equations? What is the proper partition of sets of nonclassical 
symmetry operators different from the conventional one? What are possible causes for the 
existence of singular cases for nonclassical symmetry operators? The answer to these questions 
has some fundamental importance in the research of nonclassical symmetries and will make 
finding an optimal way of obtaining the determining equation for nonclassical symmetries become 
possible. They gave a detail investigation on these questions and present a novel framework of 
singular reduction operators to clarify the main idea |57j . Here and below, following [57] we 
use sometimes the shorter and more natural term 'reduction operators' instead of 'nonclassical 
symmetry operators' or 'operators of nonclassical symmetry'. 

They also show that for any (l+l)-dimensional evolution equation [57J 

u t = H(t,x,u,ux, ...,u k ), k>l, u k = d k u/dx k , H Uk ^ 0, 

the conventional partition of the set of its reduction operators with the conditions t / and 
r = is natural since it coincides with the partition of the set into the singular and regular 
subsets. After factorizing the subsets of the reduction operator set with respect to the usual 
equivalence relation of reduction operators (see Definition [2] below) , there exist two different 
cases of inequivalent reduction operators: the regular case r = 1 and the singular case r = 
and £ = 1, which should be investigated separately. However, this is a specific property of evo- 
lution equations which does not hold for general partial differential equations in two independent 
variables. In particular, they show that for the class of nonlinear wave equations 

Utt ~ Uxx = F(u), 

where F is an arbitrary smooth function of u, which possesses two singular sets of reduction 
operators, singled out by the conditions r = £ and r = — £, and one regular set of reduction 
operators, associated with the condition r ^ ±£. The singular sets are mapped to each other 
by alternating the sign of x and hence one of them can be excluded from the consideration. 
After factorization with respect to the equivalence relation of vector fields, there are two cases 
for further study: the singular case r = £ = 1 and the regular case r ^ 1,( = 1. 

However, for more general nonlinear wave equations there exist no general results. Therefore, 
it is shown that the structures of condition symmetries of hyperbolic type nonlinear partial 
differential is more complicated than general evolution equation. 

In what follows we extend this new framework of singular reduction operators to the (1+1)- 
dimensional variable coefficient nonlinear telegraph equations ([1]). With the aid of the transfor- 
mation t = t,x = f -t\,u = u from G~/G~ 9 in section O we can reduce equation ([Tl) to one 

J 9\ x ) 

which has the same form with equation f)23f) . Thus, without loss of generality we can restrict 
ourselves to investigation in detail the equation 

f{x)u tt - {H{u)u x ) x - h{x)K{u)u x = 0, (49) 

where / = f(x),h = h(x),H = H(u) and K = K(u) are arbitrary smooth functions of the 
corresponding variables, fH > 0. 

For the sake of completeness, let us first review some necessary definitions and statements 
on nonclassical symmetries [3ll[58l[90l[96] and singular reduction operator [57] . 

6.1 Brief review of reduction operators of differential equation 

Consider an rth order differential equation C of the form L(t, x,ur r \) = for the unknown 
function u of the two independent variables t and x, where L = L[u] = L(t,x,u^ r )) is a fixed 



21 



differential function of order r and Ur r \ denotes the set of all the derivatives of the function u 
with respect to t and x of order not greater than r, including u as the derivative of order zero. 

In order to discuss the conditional symmetries of equation C, we will first treat equation C 
from a geometric point of view as an algebraic equation in the jet space J r of order r and is 
identified with the manifold of its solutions in J r [70]: C = {(t,x,Ur r \) 6 J r \L(t, x, Ur r \) = 0}. 
Let Q denote the set of vector fields of the general form 

Q = T(t,x,u)dt + £(t,x,u)d x + ri(t,x,u)d u , (t,£)^(0,0), (50) 

which is a first-order differential operator on the space M 2 x M 1 with coordinates t, x, and u. 
Then all functions invariant under Q and only such functions satisfy a first order differential 
equation 

Q[u] := ru t + £u x - r] = (51) 

called the the characteristic equation (also known as invariant surface condition). 

Denote the manifold defined by the set of all the differential consequences of the characteristic 
equation Q[u] =0 in J r by Q( r ), i-e., 

Q (r) = {{t,x,u {r) ) € J r \D?DPQ[u] = 0, a,/3eNU{0}, a + f3<r}, 

where D t = dt + u a+ i^d Ua p and D x = d x + u a ^ + id Ua p are the operators of total differentiation 
with respect to the variables t and x, and the variable u a ^ of the jet space J r corresponds to 
the derivative Q tct Q™$ ■ Denote also by Q(_ r ) * ne standard rth prolongation of Q to the space J r : 

Q(r)=Q+ ^:=rt[«]+T%fW+^l- 

0<a+^<r 

Definition 1. The differential equation C is called conditionally invariant with respect to the 
operator Q if the relation 



Q{r)[L(t,X,U(r))] 



= (52) 

£nQ (r) 



holds, which is called the conditional invariance criterion. Then Q is called conditional symmetry 
(or nonclassical symmetry, Q-conditional symmetries or reduction operator) of the equation C. 

We denote the set of reduction operators of the equation C by Q(C) which is a subset of 
Q. Any Lie symmetry operator of C belongs to Q(C). Sometimes Q(C) is exhausted by the 
operators equivalent to Lie symmetry ones in the sense of the following definition. 

Definition 2. Two differential operators Q and Q in Q are called equivalent {Q ~ Q) if 
they differ by a multiplier which is a non- vanishing function of t,x and u : Q = XQ, where 
A = \(t,x,u),\ / 0. 

Factoring Q with respect to this equivalence relation we arrive at Qf. Elements of Qf will be 
identified with their representatives in Q. The property of conditional invariance is compatible 
with this equivalence relation on Q [3T|I96| : 

Lemma 1. If the equation C is conditionally invariant with respect to the operator Q then it is 
conditionally invariant with respect to any operator which is equivalent to Q. 
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In view of this lemma, we can see that Q £ Q(C) and Q ~ Q imply Q £ Q(£), he. Q(C) 
is closed under the equivalence relation on Q. Therefore, the equivalence relation on Q induces 
a well-defined equivalence relation on Q(C); and the factorization of Q with respect to this 
equivalence relation can be naturally restricted to Q(C) that results in the subset Q/(£) of Qf. 
As in the whole set Qf, we identify elements of Qf(C) with their representatives in Q{C). In 
this approach the problem of completely describing all reduction operators for C is equivalent to 
finding Qf(C). In fact, nonclassical symmetries should be studied up to the above equivalence 
relation. The elements of Q{C) which are not equivalent to Lie invariance operators of C will 
be called pure nonclassical symmetries of C. 

The conditional invariance criterion admits the following useful reformulation |96j . 

Lemma 2. Given a differential equation C : L[u] = of order r and differential functions L[u] 
and X[u] 7^ of an order not greater than r such that L\q,< = AL|g {r) , an operator Q is a 
reduction operator of C if and only if the relation Q(r)L\cnQ( r ) = holds, where r = ordL < r 
is the order of the differential function L[u] and the manifold C is defined in J r by the equation 
L[u] = 0. 

Consider a vector field Q in the form (|5Q|) and a differential function L = L[u] of order 
ordL = r (i.e., a smooth function of variables t,x,u and derivatives of u of orders up to r). 

Definition 3. The vector field Q is called singular for the differential function L if there exists 
a differential function L = L[u] of an order less than r such that L|g, , = L\Qm ■ Otherwise Q 
is called a regular vector field for the differential function L. If the minimal order of differential 
functions whose restrictions on Q(r) coincide with L|g, , equals k (k < r) then the vector field 
Q is said to be of singularity co-order k for the differential function L. The vector field Q is 
called ultra- singular for the differential function L if L\q,, = 0. 

For convenience, the singularity co-order of ultra-singular vector fields and the order of iden- 
tically vanishing differential functions are defined to equal —1. Regular vector fields for the 
differential function L are defined to have singularity co-order r = ordL. The singularity co- 
order of a vector field Q for a differential function L will be denoted by sco^ Q. 

If Q is a singular vector field for L then any vector field equivalent to Q is singular for L with 
the same co-order of singularity. 

We will say that a vector field Q is (strongly) singular for a differential equation C if it is sin- 
gular for the differential function L[u] which is the left hand side of the canonical representation 
L[u] =0 of the equation C. Usually we will omit the attribute "strongly". 

Since left hand sides of differential equations are defined up to multipliers which are nonvan- 
ishing differential functions, the conditions from Definition [3] can be weakened when considering 
differential equations. 

Definition 4. A vector field Q is called weakly singular for the differential equation C : L[u] = 
if there exist a differential function L = L[u] of an order less than r and a nonvanishing differential 
function A = A [it] of an order not greater than r such that L\q,. = AL|g {r) . Otherwise Q is called 
a weakly regular vector field for the differential equation C. If the minimal order of differential 
functions whose restrictions on Q(r) coincide, up to nonvanishing functional multipliers, with 
L\q,. is equal to k(k < r) then the vector field Q is said to be weakly singular of co-order k for 
the differential equation C. 

The notions of ultra-singularity in the weak and the strong sense coincide. Analogous to the 
case of strong regularity, weakly regular vector fields for the differential equation C are defined 
to have weak singularity co-order r = ordL. The weak singularity co-order of a vector field Q 
for an equation L will be denoted by wsco£ Q. 
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Definition 5. A vector field Q is called a singular reduction operator of a differential equation 
C if Q is both a reduction operator of C and a weakly singular vector field of C. 

After the factorization of the reduction operator under the usual equivalence relation of re- 
duction operators in Definition [2] to singular and regular cases, the classification of reduction 
operators can be considerably enhanced and simplified by considering Lie symmetry and equiv- 
alence transformations of (classes of) equations. 

Lemma 3. Any point transformation of t, x and u induces a one-to-one mapping of Q into 
itself. Namely, the transformation g : t = T(t,x,u), x = X(t,x,u),u = U(t,x,u) generates the 
mapping g t . : Q — > Q such that the operator Q is mapped to the operator g*Q = f&i + + r)du, 
where f(t,x,u) = QT{t,x,u), £(t,x,u) = QX(t,x,u), fj(i,x,u) = QU(t,x,u). If Q' ~ Q then 
9*Q' ~ 9*Q- Therefore, the corresponding factorized mapping gf : Qf Qf also is well defined 
and bijective. 

Lemma [3] results in appearing equivalence relation between operators, which differs from 
usual one described in Definition [2j 

Definition 6. Two differential operators Q and Q in Q are called equivalent with respect to a 
group G of point transformations if there exists g G G for which the operators Q and g*Q are 
equivalent. We denote this equivalence by Q ~ Q mod G. 

The problem of finding reduction operators is more complicated than the similar problem 
for Lie symmetries because the first problem is reduced to the integration of an overdetermined 
system of nonlinear PDEs, whereas in the case of Lie symmetries one deals with a more overde- 
termined system of linear PDEs. The question occurs: could we use equivalence and gauging 
transformations in investigation of reduction operators as we do for finding Lie symmetries? 
The following statements give the positive answer. 

Lemma 4. Given any point transformation g of an equation C to an equation C, maps Q(C) 
to Q(C) bijectively. The same is true for the factorized mapping gj from Qf(C) to Qj(C). 

Corollary 1. Let G be the point symmetry group of an equation C Then the equivalence of 
operators with respect to the group G generates equivalence relations in Q(C) and in QAC). 

Consider the class C\s of equations Cq: L(t, x, ur r \,6) = parameterized with the parameter- 
functions 9 = 8(t, x, um). Here L is a fixed function of t, x, um and 9. The symbol 9 denotes the 
tuple of arbitrary (parametric) differential functions 9(t, x, u/ r \) = (0 1 (t, x, um), 9 k (t, x, «( r ))) 
running through the set S of solutions of the system S(t, x, ur r \, Or q \(t, x, Ur r \)) = 0. This system 
consists of differential equations on 9, where t,x and ur r \ play the role of independent variables 
and Qf q \ stands for the set of all the derivatives of 9 of order not greater than q. In what follows 
we call the functions 9 arbitrary elements. Denote the point transformation group preserving 
the form of the equations from C\s by G~. 

Let P denote the set of the pairs consisting of an equation Lq from C\s and an operator Q 
from Q(Cg). In view of LemmaHl the action of transformations from the equivalence group G~ 
on C\s and {Q{Cq)\9 G S} together with the pure equivalence relation of differential operators 
naturally generates an equivalence relation on P. 

Definition 7. Let 9,9' €S,Q G Q(C e ),Q' G Q{C 6 >). The pairs (C g ,Q) and {C e >,Q') are called 
G~ -equivalent if there exists g G G~ such that g transforms the equation £g to the equation/^/ , 
and Q' ~ g*Q. 

The classification of reduction operators with respect to G~ will be understood as the clas- 
sification in P with respect to this equivalence relation, a problem which can be investigated 
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similar to the usual group classification in classes of differential equations. Namely, we construct 
firstly the reduction operators that are defined for all values of 9. Then we classify, with respect 
to Cr~, the values of 6 for which the equation Cq admits additional reduction operators. 



6.2 Singular reduction operators 

Using the above notion and the procedure given by Kunzinger and Popovych in [57], we can 
obtain the following assertion. 

Proposition 1. A vector field Q = r(t,x,u)dt + £(t, x,u)d x + rj(t, x,u)d u is singular for the 
differential function L = f(x)utt — (H(u)u x ) x — h{x)K(u)u x if and only if £ 2 f(x) = t 2 H(u). 

Proof. Suppose that 0. According to the characteristic equation ru t + £,u x — r\ = 0, we can 
get 

rj £ 
ut = u x , 

T T 

Utt = {-)t ~ (-)tUx + [{-)u ~ {-)uU x \{ U x ) - (-) (-)x + {-)uU x ~ {-) X U X 

T T T T T T T T T T 

-{-)uU 2 x ~ {~)U XX \. 
T T 

Substituting the formulas of Utt from above formulaes into L, we obtain a differential function 
L = lf(x)(-) 2 ~ H(u)]u xx + f{x){(\ - £) t u x + - £) u u x }(^ - ^u x ) 

T { T T T T T T 

-(-)[(-)x + {-) u u x - {-) x u x - {-) u u 2 x ] \ - H u u 2 x - hKu x . 

T T T T T J 

According to the definition [3] of singular vector field, we have ordL < 2 if and only if f(x)(-) 2 — 
H(u) = 0. □ 

Therefore, for any /, h, H and K with fH > the differential function L = f(x)uu — 
(H(u)u x ) x — h(x)K(u)u x possesses exactly two set of singular vector fields in the reduced form, 
namely, S = {d t + y/H/ fd x + i)d u } and S* = {d t - \fW]Jd x +f)d u }, where fj = 2. Any singular 
vector field of L is equivalent to one of the above fields. Moreover, it is easy known that form 
Theorem Q] each equation of the form (|49p admits the discrete involutive transformation {x,K} 
and {x,h}. According to Corollary [H these two transformation generates a one-to-one mapping 
between S and S*. Hence it suffices, up to equivalence of vector fields (and permutation of x 
and —x ), to investigate only singular reduction operators from the set S. 

Proposition 2. For any variable coefficient nonlinear telegraph equations in the form ()49|) the 
differential function L = f(x)uu — (H(u)u x ) x — h(x)K(u)u x possesses exactly one set of singular 
vector fields in the reduced form, namely, S = {dt + y H/ fd x + fjd u }. 

Thus taking into accountant the conditional invariance criterion for an equation from class 
(|49l) and the operator dt + \J H/ fd x + rjd u , we can get 



Theorem 12. Every singular reduction operator of an equation from class (|49p is equivalent to 



Q = d t + y/H(u)/ f(x)d x + rj(t, x, u)d u , 
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where the real-valued function r](t,x,u) satisfies the determining equations 

(-h x K - 2f Vtu + \hKf x /f + - H uVVu f/H - 2f VVuu + - A fH 2 u r, 2 /H 2 

-\tlf xx /f - ^f V 2 H uu /H)^W/f- H uVx - hr)K u - ^vHufx/f - 2Hrj xu = 0, (53) 
-r]r] x H u + \/H/f(frj tt - r\ xx H - hr\ x K + 2frjrj tu + frj 2 r] uu ) = 0. 

6.3 Regular reduction operators 

The above investigation of singular reduction operators of nonlinear telegraph equation of the 
form (|49p shows that for these equation the regular case of the natural partition of the corre- 
sponding sets of reduction operators is singled out by the conditions £ ^ ±y/ 'H(u)/ fr. After 
factorization with respect to the equivalence relation of vector fields, we obtain the defining 
conditions of regular subset of reduction operator: r = 1,£ ^ ± W H (it)//- Hence we have 

Proposition 3. For any variable coefficient nonlinear telegraph equations in the form (|49p the 
differential function L = f(x)uu — (H(u)u x ) x — h{x)K(u)u x possesses exactly one set of regular 
vector fields in the reduced form, namely, S = {dt + £<9 X + fjd u } with £ ^ H(u) / f . 

Consider the conditional invariance criterion for an equation from class ()49|) and the operator 
dt + £,&x + V^u with £ =7^ iy 7 H{u)/ f , we can get 

Theorem 13. Every regular reduction operator of an equation from class ()49|) is equivalent to 

Q = d t + i{t,x)d x + r](t,x,u)d u , £(t,x) / ±y/H{u)/f 
where the real-valued functions £(t,x),r)(t,x,u) satisfy the determining equations 

UxH/f - V H U + 2H£ X + 2m = 0, 

2H U £ X + fr] uu £ 2 - r] u H u - r]H uu - Hr] uu + £H u f x /f = 0, 

-mVVuu + thKf x /f - h V K u - 2H uVx + hK&t + 2f& x - £h x K (54) 

-2ft ntu - 2fi m u ~ 2H Vxu - fiu + Hi xx = 0, 

-hKr] x - Hrj xx + fiju + 2fr] tu r] + frj 2 r] uu - 2f£ t r]x = 0. 

Solving the above system with respect to the coefficient functions £, ij, f, H, h and K under 
the equivalence group , we can get a classification of regular reduction operator for the class 
(|49p . However, due to the strong nonlinearity of system ()54p . it is difficult to get an explicit 
classification. Hence, we omit the detail investigation for the general case and concentrate on 
some special cases. 

Example. We study the regular reduction operator of equations (|49l) with H(u) = K(u) = 
u,f(x) = h(x) = 1, i.e, nonlinear telegraph equations 

u tt = (uu x ) x + uu x . (55) 

From table 3, it is easy to know that equation (|55p admits three-dimensional Lie algebra g of 
its infinitesimal Lie symmetries with a basis: 

X 1 = d t , X 2 = td t -2ud u , X 3 = d x . (56) 

The corresponding one-parameter groups are time translations and scale transformations. 
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We first discuss a special case of the regular reduction operator, i.e., consider the conditional 
symmetry operator in the form: 



Q = d x + r](t,x,u)d u . (57) 

With the assumptions r = 0, £ = 1 the determining equations (|52p for the nonlinear telegraph 
equation (jl9j) are as follows: 

Vtu = 0, r/ uu = 0, 

fVtt ~ H uuV 3 + H u {rj 2 f x /f - 3 mx - 2r?n u ) (58) 
-H(rj xx + 2r] xu r] - r)r) u f x / f - rj x f x /f) - K u hr] 2 + K(hr]f x /f - h x r] - hrj x ) = 0. 

Prom the first two equations we obtain that 

rj(t,x,u) = A(x)u + B(t,x). (59) 

Substituting the latest equation with H(u) = K(u) = u, f(x) = g(x) = 1 into the last equation 
of system (]58p . we can see that the functions A(x) and B(t, x) satisfy the overdetermined system: 



A xx + 5AA X +A X + 2A 3 + A 2 = 0, 
B xx + 5A X B + 3AB X + B X + AA 2 B + 2AB = 0, (60) 

tt 



B n - 3BB r - 2AB 2 - B 2 = 0. 



The last two equations of system (I60p imply the compatibility condition 
(UA + 4)B 2 + (36A 2 + 46A X + 18A)BB X + (2A - 37AA X - 4A 



-18.4 3 - 5A 2 )B 2 = (61) 



obtained by cross-differential. The second equation of ([60]) is a differential consequence of ([6T 
provided the equation 



(60^ - 48A 2 - 3AA - 8)B 2 - (510AA X + 122^ + 3A8A 3 + 23AA 2 + 3QA)BB X 
+(-267A 2 x - 233A 2 A X - 135AA X + 6A X - 99A 3 - 32A 2 - 70A 4 ) = 

is satisfied. 

Eliminate B x from (|6ip and (|62|) we come to the equation 



(62) 



£ 2 (272 AA X + 215 A 3 + 591 A X A 2 + 262 A 4 - 2 A + 10 A x + 230 A 2 + 38 A 2 )(1872 
-32A 2 X + 3336 A 4 - 148608 A 5 - 194152 A,.^ 3 - 11392 AA 2 X + 115668656 A 6 A x 
+363896822 A 7 A x + 11989256 A 5 A x + 539345898 A 5 A 2 X - 1263352 A 4 A x 
+5505328 A 3 A 2 X + 107208 A 3 X A + 37404240 A 3 X A 2 + 119385824 A 2 A 4 

-942408^ 2 ^ 2 + 301875108 ^ 3 ^ 3 + 4543732844* + 508439090 4 8 4 x (63) 
+979255660 4 6 42 + 800801280 4 3 4 4 + 247114530 4 2 4 4 + 638696765 4 7 42 
+695655030 4 5 4 3 + 337735335 4 4 4 3 + 59322060 44^ + 268692620 A X A 9 
+34451808 4 8 + 83063493 4 9 - 238410 4 6 + 5772399 A 7 + 6096 4 3 + 686688 A\ 
+94806420 4 10 + 16949160 4^ + 42378980 4 11 - 128 4 2 + 14960 4^42 + 128 44 x ) = 0. 



If we take the third factor in (|63p and the first equation for the function A{x) in system 
(I60p . then that overdetermined system for the function 4(x) admits solution A(x) = 0. Hence, 
solving system (I60p with A{x) = 0, we can obtain 

B(t,x) = a(t), 

where the function a(t) is given by: 
,a(t) 3 

-—==du -t-c 2 = (64) 
V6u J + oci 
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Taking a special form a(t) = 6/t 2 from (|64p yield an exact explicit solution of equation (|55p 
obtainable by solving equation (|5ip . which is an ordinary differential equation in the variable x 
for the symmetries of second type, and subsequent solution of equation (|55p for the 'constants' 
of integration actually depending on the variable t: 

150 x + 25c 2 + 25ci t 5 - 180 In (t) - 36 



u(t, x) 



25t 2 



where Ci(i = 1,2) are parameters. 

If we set B(t,x) = 0, then last two equations of (|60|) are satisfied and we arrive at an 
infinitesimal conditional symmetry 

Q = d x + A{x)ud u 

with the function A(x) satisfying the ODE A xx + 5AA X + A x + 2A 3 + ^ 2 = 0. Particular 
solution A(x) = l/4[tanh(l/2a;) — 1] of the latter equation yields exact solutions of the nonlinear 
telegraph equation ([55]) 



1 1 x 
u(t,x) = [-cosh(-x)]2 exp(--)(cii + c 2 ), 

while particular solution A(x) = l/4[coth(l/2x) — 1] yields the exact solution 

1 1 x 
u(t,x) = [sinh(-x)]5 exp(--)(cit + c 2 ). 

Finally, if we take the second factor 272 AA X + 215 A 3 + 591 A X A 2 + 262 A 4 - 2 A + 10 A x + 
230,4.2 + 38 A 2 = in ([63]), which together with the first equation of ([60]) imply A(x) = -1/2. 
Thus, we arrive at an infinitesimal conditional symmetry in the form 
u 

Q = d x + [-- + at + P]d u , 

where a,/3 are arbitrary constants. Solving equation (|5ip . we obtain an exact solutions 

u (t, x) = 2 a t + 2 p + e~ 1 l 2x c 1 AiryAi (-1/2 22/ ^ 2 "/ 3 +/3) ) 
+e- 1 / 2x c 2 AiryBi (-1/2 2VZ % +P) 

of equation (|55p . where AiryAi{t) and AiryBi{t) are the associated Airy functions of the first 
kind and the second kind respectively. 

Now we consider the conditional symmetry operator of equation (|55p in the form: 



Q = d t + £(t,x)d x +r)(t,x,u)d u . (65) 

In this way, solving the determining equations (|54p with H (u) = K{u) = u, f(x) = h{x) = 1, 
we can get £ = const, r/ = 0. 

From the above discussion we can arrival at 

Theorem 14. Equation (I49p -ff(ii) = -fT(u) = u,f(x) = h{x) = 1 is conditionally invariant 
under the following operators: 

(1) Q = d x + a(t)d u ; 

(2) Q = d x + l/4[tanh(l/2a;) - l]ud u ; 

(3) Q = d x + l/4[coth(l/2x) - l]ud u ; 

(4) Q = d x + [-% + at + P]d u ; 

(5) Q = a t + 9 X ; 

where a(t) is given by equation ([6 
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7 Conservation laws 



Apart from exact solutions, classical and nonclassical symmetry classifications, we know that 
another important subject of group analysis is the construction of conservation laws of (systems 
of) differential equations, which play an important role in mathematical physics [70J. In fact, 
the knowledge of conservation laws is useful in the numerical integration of partial differential 
equations, for example, to control numerical errors. Also, the investigation of conservation laws 
of the Korteweg-de Vries equation became a starting point of the discovery of a number of 
techniques (such as Lax pair, inverse scattering transformation, bi-Hamiltonian structures, etc.) 
to solve nonlinear evolution equations. The existence of a large number of conservation laws 
of a evolutionary partial differential equation (system) is a strong indication of its integrability. 
Conservation laws have also significant uses in the theory of non-classical transformations and 
in the theory of normal forms and asymptotic integrability. 

In this section we classify local conservation laws of equations ([I]) with characteristics depend- 
ing, at mostly, on t, x and u. For classification we use the direct method described in [49|l79]. 
To begin with, we adduce a necessary theoretical background on conservation laws, following, 
e.g., [3911701 179] and considering for simplicity the case of two independent variables t and x. See 
the above references for the general case. 

Let W be a system W(t, x, ut p \) = of I PDEs W l = 0, . . . , W l = for m unknown functions 
u = (u , . . . , u m ) of two independent variables t and x. Here ur p ) denotes the set of all the partial 
derivatives of the functions u of order not greater than p, including u as the derivatives of the 
zero order. Let W(k) denote the set of all algebraically independent differential consequences 
that have, as differential equations, orders no greater than k. We identify W/k) with the manifold 
determined by W(k) i n the jet space jw. 

Definition 8. A conserved vector of the system W is a 2-tuple F = (F (t, x, um), F 2 (t, x, um)) 
for which the divergence Divi* 1 := DtF 1 + D X F 2 vanishes for all solutions of W (i.e., Divi^^z = 
0). 

In Definition [8] and later, Dt and D x denotes the operator of total differentiation with respect 
to the variables t and x respectively. The notation V\yy means that values of V are considered 
only on solutions of the system W. 

The crucial notion of the theory of conservation laws is one of triviality and equivalence of 
conservation laws. 

Definition 9. A conserved vector F is called trivial if F l = F l + F l , i = 1, 2, where F l and F l 
are functions of t, x and derivatives of u (i.e., differential functions), F % vanish on the solutions 
of W, and the 2-tuple F = (F 1 ,F n ) is a null divergence (i.e., its divergence vanishes identically). 

Definition 10. Two conserved vectors F and F' are called equivalent if the vector-function 
F' — F is a trivial conserved vector. 

The notion of linear dependence of conserved vectors is introduced in a similar way. Namely, 
a set of conserved vectors is linearly dependent iff a linear combination of them is a trivial 
conserved vector. 

The above definitions of triviality and equivalence of conserved vectors are natural in view 
of the usual "empiric" definition of conservation laws of a system of differential equations as 
divergences of its conserved vectors, i.e. divergence expressions which vanish for all solutions 
of this system. For example, equivalent conserved vectors correspond to the same conservation 
law. However, for deeper understanding of the problem and absolutely correct calculations a 
more rigorous definition of conservation laws should be used. 

For any system W of differential equations the set CV(W) of conserved vectors of its con- 
servation laws is a linear space, and the subset CVo(W) of trivial conserved vectors is a linear 
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subspace in CV(VW). The factor space CL(W) = CV(W)/CV (W) coincides with the set of 
equivalence classes of CV(W) with respect to the equivalence relation adduced in definition [TOl 

Definition 11. The elements of CL(W) are called conservation laws of the system W, and the 
whole factor space CL(W) is called the space of conservation laws of W. 

That is why description of the set of conservation laws can be assumed as finding CL(W) 
that is equivalent to construction of either a basis if dim CL(W) < oo or a system of generatrices 
in the infinite dimensional case. The elements of CV(W) which belong to the same equivalence 
class giving a conservation law T are considered all as conserved vectors of this conservation law, 
and we will additionally identify elements from CL(W) with their representatives in CV(W). 
For F E CV(W) and T E CL(W) the notation F E T will denote that F is a conserved vector 
corresponding to the conservation law T . In contrast to the order rp of a conserved vector F as 
the maximal order of derivatives explicitly appearing in F, the order of the conservation law J- 
is called minjrp \ F E J-}. Under linear dependence of conservation laws we understand linear 
dependence of them as elements of CL(W). Therefore, in the framework of "representative" 
approach conservation laws of a system W are considered as linearly dependent if there exists 
linear combination of their representatives, which is a trivial conserved vector. 

Let the system W be totally nondegenerate [70] . Then application of the Hadamard lemma 
to the definition of conservation law and integrating by parts imply that the left hand side of 
any conservation law of W can be always presented up to the equivalence relation as a linear 
combination of left hand sides of independent equations from W with coefficients A M being 
functions of t, x and derivatives of u: 



Definition 12. Formula (166j) and the /-tuple A = (A 1 , . . . , A') are called the characteristic form 
and the characteristic of the conservation law DtF 1 + D X F 2 = correspondingly. 

The characteristic A is trivial if it vanishes for all solutions of W. Since W is nondegener- 
ate, the characteristics A and A satisfy (f66j) for the same conserved vector F and, therefore, 
are called equivalent iff A — A is a trivial characteristic. Similarly to conserved vectors, the 
set Ch(W) of characteristics corresponding to conservation laws of the system W is a linear 
space, and the subset Cho(W) of trivial characteristics is a linear subspace in Ch(W). The 
factor space Chf(W) = Ch(W)/ Cho(W) coincides with the set of equivalence classes of Ch(W) 
with respect to the above characteristic equivalence relation. 

In most of cases, we can essentially simplify and order classification of conservation laws, 
taking into account additionally symmetry transformations of a system or equivalence transfor- 
mations of a whole class of systems. Such problem is similar to one of group classification of 
differential equations. 

Proposition 4. Any point transformation g maps a class of equations in the conserved form 
into itself. More exactly, the transformation g: t = t 9 (t,x,u), x = x 9 (t,x,u), u = u 9 (t,x,u) 
prolonged to the jet space transforms the equation DtF 1 + D X F 2 = to the equation 
DtF 1 + D X F 2 = 0. The transformed conserved vector F g = (Fg,Fg) is determined by the for- 



D t F l + D X F 2 = A 1 !^ 1 + • • • + X l W l . 



(66) 



mula 



F 2 (x,u (r) ) 



F g (x,u (r) ) 



F 1 (x,U( r) )D t t + F 2 (x,u {r) )D x t 
D t t D x x - D x i D t x 

F 1 (x,u^)D t x + F 2 (x,U( r - ) )D x x 



(67) 



Dtt D x x — D x t Dtx 
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Remark 7. In the case of one dependent variable (m = 1) g can be a contact transformation: 
t = t 9 (t, x, «(!)), x = x 9 (t, x, U(i)), tt(i) = a;,U(i)). Similar note is true for the below 

statement. 

Definition 13. Let G be a symmetry group of the system W. Two conservation laws with the 
conserved vectors F and F' are called G-equivalent if there exists a transformation g E G such 
that the conserved vectors F g and are equivalent in the sense of Definition [TOl 

Any transformation g £ G induces a linear one-to-one mapping in CV(W), transforms 
trivial conserved vectors only to trivial ones (i.e. CVo(W) is invariant with respect to g* ) and 
therefore induces a linear one-to-one mapping gj in CL(W). It is obvious that 5/ preserves linear 
(in)dependence of elements in CL(W) and maps a basis (a set of generatrices) of CL(W) in a 
basis (a set of generatrices) of the same space. In such way we can consider the G-equivalence 
relation of conservation laws as well-determined on CL(W) and use it to classify conservation 
laws. 

Proposition 5. If system W admits a one-parameter group of transformations, then the in- 
finitesimal generator Q = C9i+f] a d u a of this group can be used for construction of new conserva- 
tion laws from known ones. Namely, differentiating equation (J67]) with respect to the parameter 
e and taking the value e = 0, we obtain the new conserved vector 

F = -Q {r) F l + (D^F* - {Dj?)F. (68) 

Here Qu.) denotes the r-th prolongation of the operator Q. 

Remark 8. Formula (168p can be directly extended to generalized symmetry operators (see, 
for example, [4]). A similar statement for generalized symmetry operators in evolutionary form 
(£ l = 0) was known earlier [13, 30]. It was used in [20] to introduce a notion of basis of 
conservation laws as a set which generates a whole set of conservation laws with action of 
generalized symmetry operators and operation of linear combination. 

Proposition 6. Any point transformation g between systems W and W induces a linear one- 
to-one mapping g* from CV(W) into CV(W), which maps CVo(W) into CVo(W) and generates 
a linear one-to-one mapping g^ from CL(W) into CL(VV). 

Corollary 2. Any point transformation g between systems W andW induces a linear one-to-one 
mapping g~f from Chf(W) into Chf(W). 

Consider the class W\s of systems Wq: W(t, x, ur p \ : 9) = parameterized with the parameter- 
functions = 9(t, x, Ufp-\), Here W is a fixed function of t, x, ui p \ and 0. The symbol 9 denotes the 
tuple of arbitrary (parametric) differential functions 9(t, x, U( p )) = (9 l (t, x, ut p )), ■■■,9 k (t, x, ur p \)) 
running through the set S of solutions of the system S(t, x, ur p -s, 9r q \(t, x, U( p ))) = 0. This system 
consists of differential equations on 9, where t, x and ui p \ play the role of independent variables 
and 9r q \ stands for the set of all the derivatives of 9 of order not greater than q. In what follows 
we call the functions 9 arbitrary elements. Denote the point transformation group preserving 
the form of the equations from W\s by G~ = G~(W, S). 

Consider the set P = P(W, S) of all pairs each of which consists of a system from W|s 
and a conservation law F of this system. In view of Proposition [6j action of transformations 
from G~ on W\s and {C\(We)\9 E 5} together with the pure equivalence relation of conserved 
vectors naturally generates an equivalence relation on P. 

Definition 14. Let 9,9' E 5, J 7 E CL(W g ), J 7 ' E CL(W g >),F eJ.Fe F' . The pairs (W e ,F) 
and (yVe',F') are called G~-equivalent if there exists a transformation g E G~ which transforms 
the system We to the system We' and such that the conserved vectors F g and F' are equivalent 
in the sense of Definition [TUl 
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In such a way, classification of conservation laws with respect to G7~ will be understood 
as classification in P with respect to the above equivalence relation. This problem can be 
investigated in the way that is similar to group classification in classes of systems of differential 
equations, especially it is formulated in terms of characteristics. Namely, we construct firstly the 
conservation laws that are defined for all values of the arbitrary elements. (The corresponding 
conserved vectors may depend on the arbitrary elements.) Then we classify, with respect to the 
equivalence group, arbitrary elements for each of that the system admits additional conservation 
laws. 

For more detail and rigorous proof of the correctness of the above definitions and statements 

see naiEa]. 

In what follows, we use the most direct method described in [79] to derive the conservation 
law of class (H|). Due to using the transformation t = t,x = J ^§y 5 u = u from G~ /G~ 9 in section 
El we can reduce equation ([1]) to one which has the same form with equation (I23f> . Thus, without 
loss of generality we can restrict ourselves to investigation conservation law of the equation (|23p . 

Theorem 15. A complete list of -inequivalent equations (|23p having nontrivial conservation 
laws with characteristics of the zeroth order is exhausted by ones given in table \10l 



Table 10. Conservation laws of equations (|23[) 



N 


H(u) 


K(u) 




h(x) 


Basis conservation laws 


1 


V 


V 


V 


1 


CL\ CL 2 


2 


1 






V 


CL 3 , CL 4 


3 


V 


eH + 1 


-hy 


V 


CL 5 , CL 6 


4 


V 


eH + 1 


-h y - hy' 1 


V 


CL 7 , CL 8 


5 


V 


eH + 1 


-hoj- 1 


u~ 1/2 exp ( - / |(aoo + aii)cj -1 ) 


CL 9 , CL 10 , CL 11 , CL 12 


6 


V 


eH 


V 


V 


CL 13 , CL 14 , CL 15 , CL 16 


7 


1 





V 


V 


CL 17 



Here the conserved densities F 1 and fluxes F 2 of the presented conservation laws have the 
following forms: 

CL 1 : fu t ,-(Hu x + fK); 

CL 2 : f(tu t -u), -t(Hu x + f K); 

CL 3 : e~* (h- 1 )^ + it), -e~* [hr 1 ^ - (h^u + / K}; 

CL 4 : e^hr^Uut- u), -e^/T 1 ^ - (h~%u + J K]; 

CL 5 : -e t - e f h h y (u t - u), -e t { e - e f h Hu y + hu); 

CL 6 : -e-^^hyiut + u), -e^^ 6 f h Hu y + hu); 

CL 7 : -e^^yQhi + y- l h){u t - u), -e*(e" e S h yHu y -JH + yhu); 

CL 8 : -e- t - t S h y{h y + y- x h){u t + u), -e'^e' 6 f h yHu y - f H + yhu); 

CL 9 : e- e I h f[{a n y + a w )u t - (afy + aj°)u], -{a n y + a w )(e~ e I h Hu y + hu) + a 11 J H; 

CL 10 : e- e f h f[(a 2l y + a 20 )u t - (afy + a 20 )u], -(a 2l y + a 20 )( e - e f h Hu y + hu) + a 21 J H; 

CL 11 : e- e l' h f[(a 3l y + a 30 )u t - {afy + af)u], -(a 31 y + a 30 )(e~ e f h Hu y + hu) + a 31 / H; 

CL 12 : e- e f h f[(a^y + a m )u t - (afy + af)u], -(a 41 y + a m ){e-'S h Hu y + hu) + a 41 J H; 

CL 13 : e~ e f h fu t , -e- e l' h Hu y ; 

CL 14 : e- e $ h f{tu t - it), -e" e J ' h tHu y ; 

CL 15 : e- e l' h yfu t , -e- e f h yHu y + J H; 

CL 16 : e- e l' h yf{tu t -u), -e~ e J h tyHu y + t f H; 
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CL 17 : f(au t - a t u), -au x + a x u; 

the variable y is implicitly determined by the formula x = f e~ e f h ( y ^ dy dy; e,a,ij = const, i,j = 

0,1; (a kl ,a k0 ) = (a kl (t) , a k0 (t)) , k = 1,2,3,4 is one of fundamental solution system of the 
system of ODEs 

a tt = -(an" 1 + aoi"o), oP tt = -(aioa 1 + a oao) 

w = a 0iV 2 + («00 — an)y ~~ a io! o~ = cr(t,x) is an arbitrary solution of the linear equation 
fau — o~ xx = 0. Hereafter J H = J Hdu, J K = J Kdu. In case [6] e E {0, 1} mod G±. 

In Theorem 1151 the classification of conservation laws is performed with respect to the usual 
equivalence group G± and thus some obtained results are explicit. Hence, we can further for- 
mulate the classification result in an implicit form, and indeed we can split case [5] of table [101 
into a number of inequivalent cases depending on values of . At the same time, using the 
extended equivalence group Gi , we can present the result of classification in a closed and simple 
form with a smaller number of inequivalent equations having nontrivial conservation laws. 

Theorem 16. A complete list of G± -inequivalent equations (I23p having nontrivial conservation 
laws with characteristics of the zeroth order is exhausted by ones given in table [771 



Table 11. Conservation laws of equations (|23[) 



N 


H(u) 


K(u) 


/(*) 


h(x) 


Basis conservation laws 


1 


V 


V 


V 


1 


CL 1 , CL 2 


2 


1 




h{h~ l ) xx 


V 


CL 3 , CL 4 


3 


V 


1 


-he 


V 


CL 5 , CL 6 


4 


V 


1 


— h x — hx~ x 


V 


CL 7 , CL 8 


5a 


V 





V 


V 


CL 9 , CL 10 , CL 11 , CL 12 


5b 


V 


1 


1 


1 


CL 13 , CL 14 , CL 15 , CL 16 


5c 


V 


1 


e x 


e x 


CL 17 , CL 18 , CL 19 , CL 20 


5d 


V 


1 




x» 


CL 21 , CL 22 , CL 23 , CL 24 


6a 


V 


1 


i -i 
■Are * 

X° 


ie-i 

X 


CL 25 , CL 26 , CL 27 , CL 28 


6b 


V 


1 


1 


1 

X 


CL 29 , CL 30 , CL 31 , CL 32 


7 


V 


1 


|;c — 1|^-S |ar -+- 


\x- l|"-*|a; + l|-^-i 


CL 33 , CL 34 , CL 35 , CL 36 


8 


V 


1 


g p arctan x (2 _|_ — § 


g p arctan x ^2 _|_ ]\ — 5 


CL 37 , CL 38 , CL 39 , CL 40 


9 


1 





V 


V 


CL 41 



Here the conserved densities F 1 and fluxes F 2 of the presented conservation laws have the fol- 
lowing forms: 

CL 1 : fu t , ~{Hu x + / K); 

CL 2 : f{tu t — u), —t(Hu x + J K); 

CL 3 : e-'O" 1 )^ + u), -e-^/r 1 ^ - {h~ l ) x u + / K]; 

CL 4 : e\h- l ) xx {u t - u), -e^/T 1 ^ - {h~ l ) x u + / K); 

CL 5 : e t f(u t — u), —e t (Hu x + h f K); 

CL 6 : e-*f(ut + u), -e~\Hu x + h f K); 

CL 7 : e t xf(u t — u), —e t (xHu x — J H + xhu); 

CL 8 : e~ l xf(ut + u), —e~ t (xHu x — f H + xhu); 

CL 9 : /n t , -H^; 

CL 10 : f(tu t -u), -tHu x ; 

CL 11 : ar/iij, -xHu x + J H; 
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CL 12 : xf(tu t — u), —t(xHu x — f H); 

CL 13 : ut, —Hu x — u; 

CL 14 : tut — u, —t(Hu x + u); 

CL 15 : (x - \t 2 )u t + tu, -(x- \t 2 ){Hu x + u) + J H; 

CL 16 : (tx - lt 3 )u t - (z - \t 2 )u, -{tx - \t z )(Hu x +u)+tfH; 

CL 17 : e x (ut cost + usint), — cos t(Hu x + e x u); 

CL 18 : e x (ut sint — ucost), — sint(Hu x + e x u); 

CL 19 : (x sint + cost)e x ut — (x cost + ^ cost)e x u, —(x sint + it cos t)(Hu x + e x u) + sint J ff; 

CL 20 : [(a; — |) cost— \tsm.t]e x ut + {x suit + \t cost) e x u, — [{x — \) cost— |t sint] (Hu x + e x u) + 
cost J i?; 

CL 21 : x^~ l [utCOs(^f]it) + yJJiusin{^J~(It)], — cos(^fjlt)(Hu x + x^u); 
CL 22 : x M_1 [ttt sin( v /^It) — v /^tucos( v /7 it)], — sin( v //Jt)(i/M a; + x M u); 
CL 23 : x^[u t cos(^/ / u + It) + ^ + lnsin( % /^ + It)], -cos(^//i + lt)(xHu x + x^ +1 u) 
+ cos{^fjT+lt) f H; 

CL 24 : x^[u t sm(^/JTTTt) - ^/J^TTucos(^/J^Tlt)}, - sin(^/JTTTt)(xHu x + x^ +1 u) 
+ sin(^ITTt) / H; 

CL 25 : -^e~^{ut cos t + usint), — cos t{xHu x + ^e~^u) + cos t J H; 
CL 26 : -^e~*(ut sint — ncos t), — sint(xHu x + ^e~*u) + sint J 

CL 27 : :Jre" [( — itx cos t + sin t)u t — {—\x cos t + ^tx sin t + cos t)u], — (— ^tx cost + smt)(Hu x + 
^e~*u) — ^tcost J 

CL 28 : j^e" [(|txsint + ^xcost + cost)ut — (^txcost — sint)u], — (^txsint + ^xcost + 
cos t)(Hu x + ^e~~u) + | (t sint + cost) j H; 
CL 29 : ^u t , -xHu x -u + J H; 
CL 30 : ^{tu t - u), -t(xff Ua . + U-JH); 

CL31; (5& + " "(^ + + itx) + it 2 / H; 

CL32: (g£ + ~ (|£ + £K - (^x + t)( Hu x + 1 «) + |t 3 fH; 

CL 33 : (x + l)/[sin(V2/i - lt)ut - V 2 ^ = 1 cos(V2u - lt)u], -{x + 1) sin(^Z/z - lt){Hu x + 
fat) + sin(V2/U - It) / H; 

CL 34 : (x + 1)/[cos(V2a* - lt)u t + y/2fi - lsm(y/2fj, - lt)u], -(x + 1) cos 2 - It) (Hu x + 
fat) + cos(V2^ - It) / ff; 

CL 35 : (sc - l)/[sin(V2/i + lt)n t - yJ2\i + 1 cos(V2u + lt)u], -(x - 1) sin(V2u + lt)(Hu x + 
fat) + sin(V2/x + It) / ff; 

CL 36 : (x - l)/[cos(V2/U + lt)n t + ^J2[i + 1 sin(V2/i + lt)tt], -(x - 1) cos(V2/i + lt)(Hu x + 
fat) + cos(^/2/U + It) / 

CL 37 : e-^^*(x - i)f{u t + y/-n-iu), -e-^ fZlt {x - i)(Hu x + fat) + e'^ 1 ^ 1 f H; 
CL 38 : ev 73 ^ 1 *^ - *)/(«* - -e v/= ^ =It (x - i)(Fn a . + hu) + e^"^"^ / H; 

CL 39 : e-v^ 7 ^*^ + i)/( Ui + y/-n + iu), -e'^+^^x + i)(Fu a; + fat) + e -v / =7^i J H ; 
CL 40 : e v /37r:R *(x + i)/(n t - y/-fj, + iu), -e^^^x + i)(Fn x + fat) + e V-M+»t J ij ; 
CL 41 : /(crit t - fT t u), -att^ + ^u; 
= const, a = a(t, x) is an arbitrary solution of the linear equation fott — &xx = 0. 

Proof. We search the first-order conservation laws for the equations from class (|23p in the form 
D t F(t, x, u, u t ,u x ) + D x G(t, x, u, u t ,u x ) = 0, (69) 
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where D t and D x are the operators of the total differentiation with respective to t and x corre- 
spondingly, namely, D t = dt + utd u + • • • , D x = d x + u x d u + • • • . Substituting the expression for 
u tt deduced from (j23|) into (j69|) and decompose the obtained equation with respect to u xt and 
u xx , we obtain 

F Ux + Gut = 0) 

F Ut j + G Ux =0, (70) 
Fu t + ^y- + G ^) u * + p uut + F t + G x = 0. 

Up to conserved vectors equivalence for the first equation of system (JTOj) , we can assume F Ux = 
G Ut = 0) which together with the second equation imply 

F = F 3 (t, x, u)ut + F 2 (t, x, u), G = -F 3 (t, x, u)—u x + G 1 ^, x, u). (71) 

Substituting these expression into the last equation of system (|7Q|) and splitting it with respect 
to the powers of u x and Ut, we obtain the system of PDEs for the functions F 3 ,F 2 and G 1 of 
the form 

F M 3 = 0, F 2 a +F t 3 = 0, 



p3 — + G\- jF 3 + F 3 ^H = 0, (72) 



— + G U 
F 2 + Gl=0. 
Solving first three equations of (T72j) yields 

F 3 = F 1 (t,x) 1 F 2 = -F t 1 u + F°(t,x), G l = {—) x J H- — J K + G\t,x). (73) 

Substituting the latter expression for F 2 and G 1 into the last equation of system (j72|) . we can 
know that the major role for classification is played by a differential consequence 

77I u pi 

(—) xx H-(—) x K-F t \ = 0. (74) 

Indeed, it is the unique classifying condition for this problem. In all classification cases we 
obtain the equation F® + G x = 0. Therefore, up to conserved vectors equivalence we can assume 
F° = G° = 0, and additionally F 1 7^ for conservation laws to be non-trivial. Thus, taking into 
accountant (|7ip and (|73p the conservation density and flux can be rewritten as 

(75) 



H F 1 f hF 1 f 

F = F\t,x)u t - Flu, G = -jF 1 ^ + (y) a J H- — J K. 



Equation (|74p implies that there exist no non-trivial conservation laws in the general case. Let 
us classify the special values of the parameter-functions for which equation ()23p possesses non- 
trivial conservation laws. There exist four different possibilities for values of H and K. 

1. dim(/i, H, K) = 3. It follows from flE} that = (F 1 / ' f) xx = {hF 1 / f) x = and 
therefore F 1 = a(t)f/h, (l/h)xx = 0, a# = i.e., obviously h £ {l,x _1 },a = c\ + C2t mod G±. 
Moreover, h = 1 ~ h = x^ 1 mod G^ (the corresponding transformation is x = ln|x| and 
/ = x 2 f, the other variables and parameter- functions are not changed). As a result, we obtain 
case [TJ 
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2. H € (l),i^e(l). Then = 1 mod and (hF 1 / f) x = 0,F t \ - (Fl/f) xx = 0, i.e., 
F 1 = a(t)f/h, where a^/a = A = const (otherwise we have case [I]) and so A = 1,/ = h(h ) xx 
mod (case El). 

3. He (1),K G (H,l). Then K G (1) mod and (F l /f) xx = 0,F t \ = -K(hF l /f) x , i.e. 
F 1 = (a 1 (t)a; + a (t))f and a\ t xf + a° t f = — i^(a 1 (x/i) x + a°/i I ). For if = we obtain case[5^, 
at once. Suppose K ^ 0. Xhen — 1 mod G^ and the dimension vn — dim(f,xf,h x , (xh) x ) 
can have only the values 2 and 3. 

If m = 3 then there exist constants a,ij,bi,i,j = 0,1, and a function = 9{x) such that 
(&o, &l) ¥= (0) 0)) dim(/, x/, 0) = 3 and ^ = a o/ + «oi^/ + M, (a^Os = aio/ + a,\\xf + b\0. 
Therefore, a\ t = —ana 1 — aoia ,a? t = — aioa 1 — aooa^&ia 1 + boa = 0, i.e. a 1 = c\e xt + 
C2e~ xt , a = cse St + C4e~ St where Ci(i = 1, • • • , 4), A, S = const and A, 5 ^ (otherwise, this case 
is reduced to a subcase of[l]), hence A, 6 = 1 mod G~. Depending on values (either vanishing 
or non- vanishing) of ci,C2 and 03,04 we obtain cases [3] and H] correspondingly. 

If m = 2 then ^ = aoo/ + aoixf, (xh) x = aio/ + a<il%f for some constants cnj,i,j = 0, 1. 
Therefore, a* t = —ana 1 — aoia°,a^ = — a^a 1 — aooa ,/ = —h/uj, h x /h = — (aoix + aoo)/w, 
where w = aoix 2 + (aoo — an)x — aio, i-e., a = a; -1 / 2 exp ( — / 5 (aoo + aii)w _1 ). As a results, 
we obtain four conservation laws with the conserved vectors 

(a^(t)x + a»°{t))fu t - (af(t)x + af(t))fu, -(o^(t)x + a^{t)){Hu x + hu) + a^ 1 J H), 

where (a^ 1 , a M °), \i = 1, • • • ,4 form a fundamental set of solutions of the system a\ t = —ana 1 — 
aoia°,a^ = — aioa 1 — aooa . Separate consideration of possible inequivalent values of the 
constants ajj leads to cases [5]b-[5jd and [6]|8j 

4. H,K e (1). Therefore, H= 1,K = mod G~ and F^—(F 1 /f) xx = 0. Let F 1 = a(t,x)f, 
we have ct^ — /cr^ = 0, which corresponding case[H 

□ 

The above conservation laws can be used for construction of potential systems, potential 
symmetries and potential conservation laws. We will present such analysis elsewhere. 

8 Conclusion and Remarks 

In summary, we have present an enhanced classical and nonclassical symmetries and conserva- 
tion laws analysis of the class of equations (pQ) in the framework of modern group analysis of 
differential equations. 

We have performed a complete and extended symmetry group classification of the class of 
equations (pQ) with the two "best" gauges g = 1 and g = h. The main results on classification 
are collected in tables HH3](g = 1) andHH6](g = h) where we list inequivalent cases of extensions 
with the corresponding Lie invariance algebras. The success in the classification and the clear 
presentation of the final results are relied heavily on the regular applications of four original 
tools presented in [37], i.e., the equivalence relation with respect to the extended equivalence 
group instead of the usual one, the choice of true gauges, furcate split and systematic usage of 
additional equivalences. Among them the first two kinds of techniques (the extended equivalence 
group and true gauges) are of crucial importance for obtaining a closed and explicit classification 
list. The extended equivalence group of class ([IJ) is the extension of the usual one with the non- 
trivial group of gauge equivalence transformations including transformations which are nonlocal 
in arbitrary elements. Neglecting this transformations leads to critical swelling and complication 
of both calculations and results. This can be seen from the classification results of equation ([1]) 
with the gauge g = 1 under the usual equivalence transformations (adduced in the Appendix). 
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As an application of the classification results, exact solutions of some classification models are 
given by using the method of classical Lie reduction. 

Nonclassical symmetries of equation ([1]) are discussed within the framework of singular re- 
duction operator. Determining equations related the general singular and regular reduction 
operators of (H|) with g = 1 are given. Several reduction operators of a special nonlinear tele- 
graph equation (equation ([I]) with H(u) = K(u) = u, f(x) = h(x) = 1) are present. This 
enabled to obtain some exact solutions of the corresponding equation which are invariant under 
these conditional symmetries. 

Using the most direct method, we have also investigated two classifications of local conser- 
vation laws up to equivalence relations which are generated by both usual and extended equiv- 
alence groups. Equivalence with respect to these groups and correct choice of gauge coefficients 
of equations play the major role for simple and clear formulation of the final results. 

It should be noted that all above results of class (H|) can be applied to the specific models listed 
in section [2j We do not discuss here because there are nothing but tendinous computations. 

One of the natural continuation for further investigation of different properties of class ([I]) 
is to perform group classifications of the variable gauges , i.e., when the value of g (1 or h) 
depends on values of other arbitrary elements, in a way similar to [H]. We can also make a 
further studies of non-Lie exact solutions of class (H|) by means of the classification of singular 
and regular reduction operators. Furthermore, the proposed classifications of local conservation 
laws of class ([T]) can be used to find all possible inequivalent potential systems and potential 
conservation laws (see [50] for detail) associated to the given system of differential equations. 
These problems will be investigated in subsequent publication. 

Acknowledgements 

D.j. Huang express his sincerely thanks to professor Nataliya Ivanova for stimulating discussions 
and correction of the results of this paper. This work was partially supported by the National 
Key Basic Research Project of China under Grant No. 2010CB126600, the National Natural 
Science Foundation of China under Grant No. 60873070, Shanghai Leading Academic Discipline 
Project No. B114, the Postdoctoral Science Foundation of China under Grant No. 20090450067, 
Shanghai Postdoctoral Science Foundation under Grant No. 09R21410600 and the Fundamental 
Research Funds for the Central Universities under Grant No. WM0911004. 

References 

[1] Ablowitz M J, Clarkson P A, Solitons, nonlinear evolution equations and inverse scattering, Cambridge, 
Cambridge University Press, 1991. 

[2] Abramenko A. A., Lagno V.I., Samoilenko A.M., Group classification of nonlinear evolution equations. II. 
Invariance under solvable local transformation groups, Differ. Equ. 2002, V.38, 502-509. 

[3] Akhatov I.Sh., Gazizov R.K. and Ibragimov N.Kh., Group classification of equation of nonlinear filtration 
Dokl. AN SSSR, 1987, V.293, 1033-1035. 

[4] Akhatov I.Sh., Gazizov R.K., Ibragimov N.Kh., Nonlocal symmetries. A heuristic approach, Itogi Nauki 
i Tekhniki, Current problems in mathematics. Newest results, 1989, V.34, 3-83 (Russian, translated in J. 
Soviet Math., 1991, V.55, 1401-1450). 

[5] Ames W.F., Nonlinear partial differential equations in engineering, V.l, New York, Academic, 1965, V.2, 
New York, Academic, 1972. 

[6] Ames W.F., Adams E. and Lohner R.J., Group properties of utt = [f(u)u x } x , Int. J. Non-Linear Mech., 
1981, V.16, 439-447. 

[7] Arrigo D.J., Group properties of u xx — u™u yy = f(u), Int. J. Non-Linear Mech., 1991, V.26, 619-629. 



37 



Barone A., Esposito F., Magee C.G. and Scott A.C., Theory and applications of the sine-Gordon equation, 
Riv. Nuovo Cimento, 1971, V.l, 227-267. 

Basarab-Horwath P, Lahno V.L, Zhdanov R.Z., The structure of Lie algebras and the classification problem 
for partial differential equations, Acta Applicandae Mathematicae, 2001, V.69, 43-94. 

Bluman G. and Anco S.C., Symmetry and integration methods for differential equations, Applied Mathe- 
matical Sciences, V.154, New- York, Springer- Verlag, 2002. 

Bluman G., Cheviakov A.F. and Ivanova N.M., Framework for nonlocally related PDE systems and nonlocal 
symmetries: Extension, simplification, and examples, J. Math. Phys., 2006, V.47, 113505. 

Bluman G., Cole J.D., The general similarity solution of the heat equation, J. Math. Mech., 1969, V.18, 
1025-1042. 

Bluman G., Kumei S., Symmetries and Differential Equations, Springer, New York, 1989. 

Bluman G. and Temuerchaolu, Comparing symmetries and conservation laws of nonlinear telegraph equa- 
tions, J. Math. Phys., 2005, V.46, 073513. 

Bluman G. and Temuerchaolu, Conservation laws for nonlinear telegraph equations, J. Math. Anal. Appi, 
2005, V.310, 459-476. 

Bluman G., Temuerchaolu and Sahadevan R., Local and nonlocal symmetries for nonlinear telegraph equa- 
tion, J. Math. Phys., 2005, V.46, 023505. 

Boyko V.M. and Popovych V.O., Group classification of Galilei-invariant higher-orders equations, Proceed- 
ings of Institute of Mathematics of NAS of Ukraine, 2001, V.36, 45-50. 

Cattaneo C, Sulla conduzione del calore, Atti Sem. Mat. Fts. Univ. Modena, 1948, V.3, 83-101. 

Chikwendu S.C., Non-linear wave propagation solutions by Fourier transform perturbation, Int. J. Non- 
Linear Mech., 1981, V.16, 117-128. 

Clarkson P. A. and Mansfield E.L., Symmetry reductions and exact solutions of a class of nonlinear heat 
equations, Physica D, 1994, V.70, 250-288. 

Clarkson P. A. and Winternitz P., Nonclassical symmetry reductions for the Kadomtsev-Petviashvili equa- 
tion, Physica D, 1991, V.49 257-272. 

Cristescu N., Dynamic Plasticity, Northholland, Amsterdam (1967). 

Donato A., Similarity analysis and nonlinear wave propagation, Int. J. Non-Linear Mech., 1987, V.22, 
307-314. 

Donato A. and Fusco D., Wave features and infinitesimal group analysis for a second order quasilinear 
equation in conservative form, Int. J. Non-Linear Mech., 1987, V.22, 37-46. 

Engelbrecht J., Nonlinear wave processes of deformation in solid, Pitman Publishing (1983). 

Faddeev L.D., Takhtajan L.A., Hamiltonian method in the theory of solitons, Berlin, Springer- Verlag, 1987. 

Foursov M.V. and Vorob'ev E.M., Solutions of the nonlinear wave equation utt = (uu x ) x invariant under 
conditional symmetries, J. Phys. A: Math. Gen., 1996, V.29, 6363-6373. 

Fushchych W.I., Shtelen W.M. and Serov N.L, Symmetry Analysis and Exact Solutions of Nonlinear Equa- 
tions of Mathematical Physics, Dordrecht: Kluwer, (English transl.) 1993. 

Fushchych W. L, Shtelen W.M., Serov M.I. and Popovych R.O., Q-conditional symmetry of the linear heat 
equation, Proc. Acad. Sci. Ukraine, 1992, no. 12, 28C33. 

Fushchych W.I., Tsyfra I.M., On a reduction and solutions of the nonlinear wave equations with broken 
symmetry, J. Phys. A: Math. Gen., 1987, V.20, L45-L48. 

Fushchych W.I. and Zhdanov R.Z., Conditional symmetry and reduction of partial differential equations, 
Ukr. Math. J., 1992, V.44, 970C982. 



38 



[32] Gandarias M.L., Torrisi M., and Valenti A., Symmetry classification and optimal systems of a non-linear 
wave equation, Int. J. Non-Linear Mech., 2004, V.39, 389-398. 

[33] Gagnon L., Winternitz P. Symmetry classes of variable coefficient nonlinear Schroinger equations. J. Phys. 
A: Math. Gen., 1993, V.26, 7061-7076. 

[34] Gazeau J P, Winternitza P. Symmetries of variable coefficient Korteweg-de Vries equations . J. Math. Phys., 
1992, V.33, 4087-4102. 

[35] Goenner H., Havas P., Exact solutions of the generalized LaneCEmden equation, J. Math. Phys., 2000, V.41, 
7029-7042. 

[36] Grundland A.M. and Tafel J., On the existence of nonclassical symmetries of partial differential equations, 
J. Math. Phys., 1995, V.36, 1426-1434. 

[37] Huang D.J., Ivanova N M, Group analysis and exact solutions of a class of variable coefficient nonlinear 
telegraph equations, J. Math. Phys., 2007, V.48, 073507. (23 pages) 

[38] Huang D.J., Mei J.Q., Zhang H.Q., Group classification and exact solutions of a class of variable coefficient 
nonlinear wave equations, Chin. Phys. Lett., 2009, V.26, 050202. 

[39] Huang D.J., Zhou S.G., Group properties of generalized quasi-linear wave equations, J. Math. Anal. Appl., 
2010, V.366, 460-472. 

[40] Huang D.J., and Zhang H.Q, Preliminary group classification of quasilinear third-order evolution equations, 
Appl. Math. Mech. 2009, V.30(3), 275-292. 

[41] Huang Q., Lahno V., Qu C.Z. and Zhdanov R., Preliminary group classification of a class of fourth-order 
evolution equations, J. Math. Phys., 2009, V.50 , 023503. 

[42] Huang Q., Qu C.Z., Zhdanov R., Nonlocal symmetries of fourth-order nonlinear evolution equations. 

\arXiv:0905.2033\ 

[43] Ibragimov N.H., Transformation groups applied to mathematical physics, Mathematics and its Applications 
(Soviet Series), Dordrecht, D. Reidel Publishing Co., 1985. 

[44] Ibragimov N.H. (Editor), Lie group analysis of differential equations — symmetries, exact solutions and 
conservation laws, V.l, Boca Raton, FL, CRC Press, 1994. 

[45] Ibragimov N.H., Elementary Lie group analysis and ordinary differential equations, New York, Wiley, 1999. 

[46] Ibragimov N.H., Torrisi M. and Valenti A., Preliminary group classification of equations v t t = f(x, v x )v xx + 
g(x,v x ), J. Math. Phys., 1991, V.32, 2988-2995. 

[47] Ivanova N.M., Popovych R.O. and Sophocleous C, Group analysis of variable coefficient diffusion-convection 
equations. I. Enhanced group classification, Lobachevskii Journal of Mathematics, 2010 V.31(2), 100-122. 
(arXiv:0710.273, [math-ph]) 

[48] Ivanova N.M., Popovych R.O. and Sophocleous C, Group analysis of variable coefficient diffusion-convection 
equations. II. Contractions and Exact Solutions. arXiv:0710.3049, [math-ph]. 

[49] Ivanova N.M., Popovych R.O. and Sophocleous C, Group analysis of variable coefficient diffusion-convection 
equations. III. Conservation Laws. arXiv.0110.3053 , [math-ph]. 

[50] Ivanova N.M., Popovych R.O. and Sophocleous C, Group analysis of variable coefficient diffusion-convection 
equations. IV. Potential Symmetries. arXiv.071 0.425 1 , [math-ph]. 

[51] Ivanova N.M. and Sophocleous C, On the group classification of variable coefficient nonlinear diffusion- 
convection equations, J. Comp. and Appl. Math., 2006, V.197, 322-344. 

[52] Kingston J.G. and Sophocleous C, On form-preserving point transformations of partial differential equa- 
tions, J. Phys. A: Math. Gen., 1998, V.31, 1597-1619. 

[53] Katayev I.G., Electromagnetic Shock Waves slliffe, London, 1966. 

[54] Keller J.B. and Ting Lu, Periodic vibations of systems governed by non-linear partial differential equations, 
Comm. Pure Appl. Math., 1966, 19(4), 371-420. 



39 



Kingston J.G. and Sophocleous C, Symmetries and form-preserving transformations of one-dimensional 
wave equations with dissipation, Int. J. Non-Lin. Mech., 2001, V.36, 987-997. 

Kumei S., Invariance transformations, invariance group transformations and invariance groups of the sine- 
Gordon equations, J. Math. Phys., 1975, V.16, 2461-2468. 

Kunzinger M. and Popovych R.O., Singular reduction operators in two dimensions, J. Phys. A, 2008, V.41, 
505201, 24 pp.. larXiv:0808.3577l 

Kunzinger M. and Popovych R.O., Is a nonclassical symmetry a symmetry, Proceedings of 4th Workshop 
Group Analysis of Differential Equations and Integrability, 2009. (math-ph/0903.0821) 

Lahno V.I. and Zhdanov R.Z., Group classification of nonlinear wave equations, J. Math. Phys., 2005, V.46, 
053301. 

Lagno V.I., Samoilenko A.M., Group classification of nonlinear evolution equations. I. Invariance under 
semisimple local transformation groups, Differ. Equ. 2002, V.38, 384-391. 

Levi D. and Winternitz P., Non-classical symmetry reduction: example of the Boussinesq equation, J. Phys. 
A: Math. Gen. 22 (1989) 2915-2924. 

Lie S., On integration of a Class of Linear Partial Differential Equations by Means of Definite Integrals, CRC 
Handbook of Lie Group Analysis of Differential Equations, V.2, 473-508. (Translation by N.H. Ibragimov of 
Arch, for Math., Bd. VI, Heft 3, 328-368, Kristiania 1881). 

Lisle I.G., Equivalence transformations for classes of differential equations, Thesis, University of British 



Columbia, 1992 (http://www.ise.canberra.edu.au/mathstat/StafTPages/LisleDissertation.pdf). (See also 



Lisle I.G. and Reid G.J., Symmetry classification using invariant moving frames, ORCCA Technical Report 
TR-00-08 (University of Western Ontario), [http://www.orcca.on.ca/TechReports/2000/TR-00-08.htmll ) 

Meleshko S.V., Group classification of equations of two-dimensional gas motions, Prikl. Mat. Mekh., 1994, 
V.58, 56-62 (in Russian) (translation in J. Appl. Math. Mech., 58 (1994) 629-635). 

Nikitin A.G. and Popovych R.O., Group classification of nonlinear Schrodinger equations, Ukr. Math. J., 
2001 V.53, 1053-1060. 

Novikov S.P., Manakov S.V., Pitaevskii L.P., Zacharov V.E., Theory of Solitons, The Inverse Scattering 
Method, New York, Consultants Bureau, 1980. 

Nucci M.C., Nonclassical symmetries as special solutions of heir-equations, J. Math. Anal. Appl, 2003, 
V.279, 168-179. 

Nucci M.C. and Clarkson P.A., The nonclassical method is more general than the direct method for symmetry 
reductions: an example of the Fitzhugh-Nagumo equation, Phys. Lett. A, 1992, V.164, 49-56. 

Nucci M.C. and Leach P.G.L., The determination of nonlocal symmetries by the technique of reduction of 
order, J. Math. Anal. Appl., 2000, V.251, 871-884. 

Olver P.J., Application of Lie Groups to Differential Equations, New York, Springer- Verlag, 1986. 

Oron A. and Rosenau P., Some symmetries of the nonlinear heat and wave equations, Phys. Lett. A, 1986, 
V.118, 172-176. 

Ovsyannikov L.V., Dokl. Akad. Nauk SSSR, 1959, V.125, 592-595. 

Ovsiannikov L.V., Group analysis of differential equations, 1982, New York: Academic Press. 

Patera J. and Winternitz P., Subalgebras of real three- and four-dimensional Lie algebras J. Math. Phys., 
1977 ,V.18, 1449-1455. 

Polyanin A.D. and Zaitsev V.F., Handbook of exact solutions for ordinary differential equations, Chapman 
Hall/CRC, Boca Raton, 2003. 

Popovych R.O., Classification of admissible transformations of differential equations, Collection of Works 
of Institute of Mathematics, Kyiv, 2006, V.3, N 2, 239-254. 



40 



[77] Popovych R.O. and Cherniha R.M., Complete classification of Lie symmetries of systems of two-dimensional 
Laplace equations, Proceedings of Institute of Mathematics of NAS of Ukraine, 2001 V.36 212-221. 

[78] Popovych R.O. and Ivanova N.M., New results on group classification of nonlinear diffusion-convection 
equations, J. Phys. A: Math. Gen., 2004, V.37, 7547-7565 (|math-ph/0306035|. 



[79] Popovych R.O. and Ivanova N.M., Hierarchy of conservation laws of diffusion-convection equations, J. Math. 
Phys., 2005, V.46, 043502 ( |math-ph/0407008| . 

[80] Popovych R.O., Ivanova N.M. and Eshraghi H. Group classification of (l+l)-dimensional Schrodinger equa- 
tions with potentials and power nonlinearities, J. Math. Phys., 2004, V.45, 3049-3057 (math-ph/0311039). 

[81] Popovych R.O., Kunzinger M. and Eshraghi H., Admissible point transformations and normalized classes 
of nonlinear Schr?dinger equations, Acta Appl. Math., 2010, V.109, 315-359, |arXiv:math-ph/061106l] 

[82] Pucci E., Group analysis of the equation utt + Xu xx = g(u,u x ), Riv. Mat. Univ. Parma, 1987, V.12 N 4, 
71-87. 

[83] Pucci E. and Salvatori M.C., Group properties of a class of semilinear hyperbolic equations, Int. J. Non- 
linear Mech., 1986, V.21, 147-155. 

[84] Rogers C, Ruggeri T., A reciprocal Backlund transformation: aplication to nonlinear hyperbolic boundary 
value problems, Lett. Nuovo Cimento, 1985, V.44, 289-296. 

[85] Stephani H., Differential equation: their solution using symmetries, Cambridge, Cambridge University Press, 
1994. 

[86] Suhubi E.S. and Bakkaloglu A., Group properties and similarity solutions for a quasi-linear wave equation 
in the plane, Int. J. Non-Linear Mech., 1991, V.26, 567-584. 

[87] Torrisi M. and Valenti A., Group properties and invariant solutions for infinitesimal transformations of a 
nonlinear wave equation, Int. J. Non-Linear Mech., 1985, V.20, 135-144. 

[88] Torrisi M. and Valenti A., Group analysis and some solutions of a nonliear wave equation, Atti Sem. Mat. 
Fis. Univ. Modena, 1990, V.XXXVIII, 445-458. 

[89] Vaneeva O.O., Johnpillai A.G., Popovych R.O. and Sophocleous C, Enhanced group analysis and conserva- 
tion laws of variable coefficient reaction-diffusion equations with power nonlinearities, J. Math. Anal. Appl. , 
2007, V.330, 1363-1386, [arXiv:math-ph/0605081 1 

[90] Vaneeva O.O., Popovych R.O. and Sophocleous C, Enhanced group analysis and exact solutions of vari- 
able coefficient semilinear diffusion equations with a power source, Acta Appl. Math., 2009, V.106, 1-46, 
larXiv:0708. 34571 

[91] Varley E., Seymour B., Exact solutions for large amplitude waves in dispersive and dissipative systems, 
Studies in Appl. Math., 1985, V.72, 241-262. 

[92] Vasilenko O.F. and Yehorchenko LA., Group classification of multidimensional nonlinear wave equations 
Proceedings of Institute of Mathematics of NAS of Ukraine, 2001, V.36, 63-66. 

[93] Zhdanov R.Z., Lahno V.I., Group classification of heat conductivity equations with a nonlinear source, 
J. Phys. A. : Math. Gen., 1999, V.32, 7405-7418. 

[94] Zhdanov R.Z. and Lahno V.I., Group classification of the general evolution equation: Local and quasilocal 
symmetries, SIGMA, 2005, V.l, 009. 

[95] Zhdanov R.Z. and Lahno V.L, Group classification of the general second-order evolution equation: semisim- 
ple invariance groups, J. Phys. A: Math. Theor., 2007, V.40, 5083-5103. 

[96] Zhdanov R.Z., Tsyfra I.M., Popovych R.O., A precise definition of reduction of partial differential equations, 
J. Math. Anal. Appl, 1999, V.238, 101-123. ( |arXiv:math-ph/020702"3"j ) 



41 



A Note on classification with respect to the usual equivalence 
group 

In tables [T2HT41 we list all possible G~-inequivalent sets of functions f(x), g(x), H(u), K{u) and 
corresponding invar iance algebras. 



Table 12. Case of Vff(u) 



N 


K{u) 


/(*) 


9{x) 


1 


V 


V 


V 


2 


V 


M A 


x- 1 


3 


V 


i 


(a + Px)- 1 


4 


H + /3 


(a + /3x)~% 


(a + Px)- 1 


5 


H 


V 


-f/2f + ay/\f\ 


6 


H 


V 


t 

2/ 


7 


H 


V 


8 


1 


1 


x- 1 



Basis of A n 



d t 

d t , \{\ + 2)td t + xd x 

d t , d x , td t + xd x 
d t ,d x ,{p-l)td t + pxd x 
dt,\f\-^d x 



dt,td t + xd x 



Here A / mod G~, a, /3 ^ 0. 



Table 13. Case of H(it) = 



/(*) 



Basis of A n 



7 

8 
9 

10 



1 
1 

e vu + hie v 



e vu +hie n 

ue u + h t e u 
e u 



- (iL 



= ( 



y_i£i _ j_ 

/ ^/\T\dx 2 / 

/'/ vlTj^y 



2/VI/l 

^=-2(g- 1 ) x g 
V 

/ 2 (*) 



1 r /'Iv / l7T<fa v 
hi t. 



2/VI/l 



V 



dt, \{jp — pv — 2v — q + l)tdt 

+ (1 - v)xd x + (q+l)d u 

dt, \(p- pv - q)td t + (1 - f)3a; + qd u 

dt, + p - q)tdt + xdx + (q + l)d u 

dt,\{p- q)td t + dx + qd u 

V I J 



dt, 1(1 - 2i/)tft + (1 - v ) s ^^ x d x + d u 



dt,-±/3td t +g- 1 d x + pd u 

d t , \td t - d u 

t = ci + §c 2 £,£ = Ci x ), 

^/ = -c2 + (-^-|3 1 )e 1 i 

t = ci + ±e 2 f,£ = £ 2 (z), 



r ? = - C2 + 2 C3 + (-4-| 5 2 )C 2 ( a; ) 



Here e = ±1 and q,a, f3 7^ 0. f 1 (x), g 1 {x) and ^(x) satisfy the relation 



$ M + (2$ - + * M + (# - 2*)*, - + *) = 0, £ - $^ = 0, 



where 

$ = 4 (2 7T +ffl) ' * = (^t-p 1 ); 

and f 2 (x),g 2 (x) and £ 2 (x) satisfy the relation 



+ a, - 9A - a 2 = 0, £-ez 2 = c :1 , 



where 



& = -\{^+g 2 ), A=(Jp-g 2 ) 
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The operators from tables [12H14I form bases of the maximal invariance algebras if the cor- 
responding sets of the functions /, g, H, K are G~-inequivalent to ones with most extensive 
invariance algebras. For example, in case 114111 (p, v) 7^ (0,0). Similarly, in case 113111 the con- 
straint set on the parameters fx, v and A coincides with the one for case 11411^ and /i = 1 if 
v = 0. 



Table 14. Case of H(u) = u" 



N 


A* 


K(u) 


/(*) 


g(x) 


Basis of A max 


1 


V 


u v 


\x\ p 


\x\" 


d t , $[(p-q + l)ii-(p + 2)v]tdt 
+ (m — u)xd x + (q + l)ud u 


2 


V 


u v 


e px 


ee qx 


dt, \ [{p - q)p - pv\d t + (p - v )d x + qud u 


3 


V 


u v + hm* 1 




( V 


dt, \p.td t + ^TTjjy^ ~ u ® u 








- c - f ' v 

2 f a^TT 














4 


V 






1 //'.f^iTl^y 


8 t , \{p - 2u)td t + (ja- v )LM^l dx + udu 






y\f\ 








-(f ) v\f\ dx y 
















5 


V 


1 




|sb|« 


dt, |/Lt(l + p — q)tdt + pxd x + (1 + q)ud u 


6 


V 


1 






d t , §m(p — q)tdt + nd x + qudu 


7 


V 


u M Inu + hm* 1 


4 = -2(fl -1 )'fl 


(•T 1 )" = 


dt, —\fi\itdt + g~ 1 d x + fiudu 


8 


V 


it" 


V 


V 


dt, i ixtdt — udu 


9 


' 3 


u" 


/ 3 W 




T = \nc 2 t + Ci, £ = f 3 (x) 
r) = [(4$ 3 + *3)C 3 - c 2 ]u 


10 


-4 

r 3 


it" 


/ 4 (z) 




T = \pCit + Cl, £ = £ 4 (x) 

?7 = [(4$ 4 + *4)C 4 - C2 + 2c 3 ]M 


11 


4 
3 


4 

U 3 


/ s (z) 


g 5 (x) 


r = c 2 t + ci,£ = c 3 /$ 5 

ry = jfj^ca^ + (c 2 - c 3 )$i + 2c 3 $53 5 ]^ 


12 


4 
3 


4 

tt 3 


V 


2/ 


r = c 2 f + ci,£ = C 6 (a:) 


13 





U 


\x\ p 




9 t) |(p + l)td t + xdx -(q+ l)ud u 


14 





U 


e px 




dt, \ (p + 2)t9 t + d x - qudu 


15 





In w + ho 


g 2 exp( / gdx) 


Or 1 )" = 


dt^tdt+g^dx+Pudu 


16 





u 


g 2 exp( / gcfa) 




dt,\td t +g- x dx+Pdu 


17 





e u + h 




(g^ 1 J gdx)" = hog 


d t ,-\{p + 2)td t - cT 1 / gdxdx + du 


18 





u" + ho 




(g' 1 J gdx)" = h g 


d t ,—\v{p + 2)td t - vg~ x J gdxdx + ud u 



Here e = ±1 and q, a, f3, ho, hi 7^ 0. f 3 (x),g 3 (x) and £ 3 (x) satisfy the relation 

■ip$ 3 xx - [(16/i + 32)$ 3 + (7/i + 8)*3]$3* + P^3xx - [(5/i + 8)$ 3 + (2/i + 2)*3]*3z + (24^ + 32)3-1 
+(14/i + 16)$|*3 + (2 M + 2)$ 3 *| =0, £ - = 0, 

where 

* = M 3 2( M + A_t 3 , (5^ + 6)/ x 3 

3 3/i + 4 3 (3/X + 4)/ 3 ' 3 3/i + 4 ff + (3 M + 4)/3 ' 

f 4 {x), g 4 (x) and £ 4 (a:) satisfy the relation 

4/i$ 4i r + M*4* - (12m + 16)^| - (7m + 8)$ 4 *4 - (n + 1)$1 = 0, ^ - m^C 4 = c 3 , 

where 

» = M 4 2(m + l)/, 4 M 4 (5M + 6)/ 4 

4 3m + 4 3 (3m + 4)/ 4 ' 4 3m + 4 ff + (3m + 4)/4 ' 

f (x),g (x) satisfy the relation 
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where 

'!> . f I -V: 

and £, e '(x) satisfies third order ordinary differential equation 

£,xxx — K^j) _ 2(— + [( — + = 0. 
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